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20.  Abstract  (continued) 


"taatches  InP  and  ZnGeAs2  which  matches  GaAs. 

Using  liquid  phase  epitaxial  methods  it  was  found  that  the 
lattice  match  between  CdSnP2  and  InP  is  not  close  enough  to  obtain 
high-quality  layers.  The  layer  morphology  was  greatly  improved  by 
growing  lattice  matched  ZnxCdi-xSnP2  on  InP.  One  of  these  layers 
had  a  room  temperature  electron  mobility  of  2000  cm2/V  sec  at  a 
concentration  of  3xl0^-^cm“2,  which  is  substantially  larger  than  can 
be  obtained  in  InP  at  the  same  concentration.  At  the  other  end  of 
this  chalcopyrite  alloy,  ZnSnP2  was  found  to  closely  lattice  match 
GaAs  with  good  morphology.  Some  of  these  layers  exhibited  hole 
mobilities  of  70  cm2/v  sec  at  concentrations  of  about  lxlO^-^cm”^. 
Attempts  to  obtain  n-type  conductivity  in  ZnSnP2  were  unsuccessful. 

For  both  CdSnP2  and  ZnSnP2  the  conductivity  increased  with  increased 
loss  of  volatile  constituents. 

Several  chloride-transport  vapor-phase  methods  were  examined 
in  an  unsuccessful  effort  to  grow  ZnGeAs2  layers  epitaxially  on 
GaAs  substrates.  Our  experience  with  these  methods  indicates  that 
most  of  the  reproducibility  problems  could  be  overcome  with  a 
hydride  transport  process. 

v  Discussions  in  the  literature  have  pointed  out  several  problems 
with  the  usual  theory  of  heterostructures.  We  present  a  general, 
first-order,  nonequilibrium  thermodynamic  analysis  which  points  out 
that  one  problem  with  conventional  theory  is  that  it  neglects  a 
generalized  thermodynamic  force.  This  force  can  be  associated  with 
the  gradient  of  the  effective  density-of-states .  In  heterostructures 
such  gradients  cause  the  charge  carriers  to  move  into  regions  with 
higher  density-of-states,  thus  increasing  the  entropy  of  the  system 
in  agreement  with  the  second  law  of  thermodynamics .  We  also  show 
that  if  the  electrochemical  potential  is  referenced  to  the  infinite 
vacuum  level,  then  the  electron  affinity  rule  must  be  valid  in 
the  absence  of  interface  states. ^ - 

One  application  for  some  of  the  II-IV-V2  chalcopyrites  may  be 
in  inertial  transport  devices.  Calculations  are  presented  on  the 
I-V  characteristics  of  a  space-charge  limited  diode.  If  the  electrons 
exhibit  inertial  motion  (are  not  scattered) ,  then  they  will  quickly 
reach  nonparabolic  regions  of  the  conduction  band.  When  this 
nonparabolicity  is  taken  into  account,  the  current  appears  ohmic 
rather  than  exhibiting  a  vV2  behavior,  as  has  been  previously 
supposed.  Inertial  transport  may  also  be  utilized  to  achieve  a 
terahertz  Zener  oscillator.  If  Zener  oscillations  are  to  be  used 
as  coherent  radiation  sources,  the  carriers  must  oscillate  in  phase. 

A  scheme  is  proposed  for  phase  initialization  and  focusing,  whereby 
carriers  are  injected  into  an  interaction  region  in  phase,  and 
scattered  carriers  are  eliminated  from  the  interaction  region.  The 
method  utilizes  injection  across  a  heterojunction  energy  barrier, 
a  submicron  interaction  region,  and  electromagnetic  feedback. 
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1 .  INTRODUCTION 

The  II-IV-V2  chalcopyrites  have  several  potential  advan¬ 
tages  over  III-V  compounds  and  alloys  for  high  speed  device 
applications.  For  example,  ZnGeAs2  and  CdSnP2  are  expected ^ 
to  have  lower  conduction  band  effective  masses  than  their 

technologically  important  III-V  compound  analogs,  GaAs  and  InP. 

.  f  21 

This  together  with  the  differences  in  band  structure  J  could 

result  in  higher  hot  electron  mobilities  and  saturated  drift 
velocities  for  the  II-IV-V2  chalcopyrites.  One  might  also  ex¬ 
pect  some  of  the  ordered  II-XV-V2  compounds  to  have  higher  carrier 
mobilities  and  thermal  conductivities  than  their  equivalent 
disordered  III-V  alloys.  The  technology  of  the  II-XV-V2  chal¬ 
copyrites,  however,  is  not  as  advanced  as  that  of  their  III-V 
analogs  and  few  of  these  potential  advantages  have  been  realized. 
In  fact,  the  state  of  the  art  for  these  materials  is  about 
the  same  as  that  for  the  III-V  compounds  a  number  of  years 
ago,  when  the  current  interest  in  III-V  compounds  for  device 
applications  began  to  develop.  Considering  the  potential  ad- 
vantages_  of  these  materials  this  work  was  undertaken  to  develop 
some  of  the  II-IV-V2  compounds  into  technologically  useful 
materials . 

1.1  MATERIAL  CONSIDERATIONS 

Just  as  it  was  for  the  III-V  compounds,  it  is  difficult 
to  determine  which  of  the  II-IV-V2  chalcopyrites  are  most 
promising  for  high-speed  device  applications.  If  we  take  into 
account  the  evolution  of  the  III-V  compounds ,  we  would  choose 
the  ternary  analogs  of  GaAs  and  InP:  ZnGeAs2  and  CdSnP2. 


r  2 1 

ZnGeAs2  has  been  determined  to  have  an  energy  gap  of  1.15  eVl  , 
somewhat  lower  than  the  1.43  eV  of  GaAs.  With  a  lattice  con¬ 
stant  of  5.672  A[3] 

it  should  be  possible  to  grow  it  epitaxially 
on  Ge  with  a  lattice  mismatch  of  0.25  percent  or  on  GaAs  with 
a  mismatch  of  0.34  percent.  However,  it  has  only  been  grown 

f  31 

from  the  melt  with  high  p-type  conductivity  and  low  mobility  . 

CdSnP2  has  been  determined  to  have  an  energy  gap  of  about 
r  4 1 

1.17  eV  ,  somewhat  lower  than  the  1.35  eV  of  InP.  With  a 

o  r  3 1 

lattice  constant  of  5.900  A  it  could  be  grown  epitaxially  on 

InP  with  a  mismatch  of  0.53  percent.  More  work  has  been  done 

on,  CdSnP2  than  on  ZnGeAs2  ^  .  CdSnP2  has  previously  been  grown 

epitaxially  on  InP  from  a  Sn  melt  with  low  n-type  conductivity 

and  reasonable  mobility  It  has  also  been  made  p-type 

The  conduction  band  effective  mass  for  CdSnP2  has  been  measured 

to  be  0.036  which  is  better  than  half  the  mass  of  InP  and 

less  than  the  calculated  value  of  0.059  m^  for  CdSnP2. 

To  utilize  the  potential  high-speed  capabilities  of  any 

of  the  II-IV-Vp  chalcopyrites  it  is  necessary  to  grow  well- 

controlled  epitaxial  layers  of  the  material  on  semi-insulating 

and  heavily-doped  substrates.  An  examination  of  the  CdP2~Sn 
r  gi 

phase  diagram  indicates  that  CdSnP2  forms  peritectically 
and  that  it  can  be  grown  from  a  Sn-rich  solution  at  temperatures 

O 

between  220  and  570  C.  As  indicated,  this  method  has  been  used 
to  grow  layers  of  CdSnP2  epitaxially  on  InP  substrates^  and 
has  also  been  used  to  grow  platelets  of  CdSnP2<  Thus,  liquid 
epitaxy  seems  to  be  a  good  technique  for  growing  CdSnP2. 


The  ZnGeAs2  phase  diagram  [9]  indicates  a  high  solubility  of  Ge 
which  precludes  the  use  of  liquid  epitaxy  to  grow  this  material. 
This  high  Ge  solubility  may  account  for  the  poor  quality  of 
material  grown  from  the  melt.  For  this  reason  a  vapor-phase 
technique  should  be  a  suitable  method  for  growing  ZnGeAs2 
epitaxially  on  GaAs  or  Ge  substrates. 

The  heterostructure  properties  of  CdSnP2  grown  on  In? 
substrates  and  ZnGeAs2  on  GaAs  or  Ge  could  be  helpful  for  high 
speed  device  applications.  The  somewhat  lower  band-gap  of 
CdSnP2  compared  to  InP  and  ZnGeAs2  compared  to  GaAs  could  help 
confine  carriers  to  the  epitaxial  layer  and  thus  prevent  the 
injection  of  hot  electrons  into  the  substrate  in  the  same 
manner  as  has  been  discussed  [10]  for  InGaAs  FETs  on  GaAs. 

Also,  to  obtain  exact  lattice  matched  heterostructures,  it  may 
be  possible  to  grow  alloys  of  these  materials  with  other  chal- 
copyrites  on  InP  and  GaAs. 

1.2  RESEARCH  OBJECTIVES 

The  objectives  of  this  work  as  they  evolved  from  these 
arguments  and  the  work  itself  were  to: 

(1)  develop  techniques  for  the  growth  of  high-quality 
epitaxial  Zn  Cd.  SnP9  on  InP  and  GaAs  and  ZnGeAs, 

X  1"X  fc  ^ 

on  GaAs  or  Ge; 

(2)  minimize  electrically-active  defects  and  residual 
impurities  in  these  materials  to  achieve  well-controlled 
doping ; 


(3)  investigate  the  electrical  transport  properties  of 
these  materials  to  determine  their  suitability  for 


heterostructures  and  high-speed  device  applications; 

(4)  develop  the  technology  required  for  the  fabrication 
of  suitable  device  structures;  and 

(5)  investigate  novel  high-speed  device  concepts  which 
could  utilize  the  unique  properties  of  these  chalco- 
pyrite  materials. 
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2.  LIQUID  PHASE  EPITAXIAL  GROWTH  OF  Zn  Cd,  SnP- 

_ _  X  1-X  4 


CdSnP2  is  a  semiconductor  material  of  some  interest 
for  high  speed  and  microwave  devices  because  it  has  a 
smaller  electron  effective  mass  than  its  III-V  analog 
InP  [1]  and  because  of  the  observation  of  negative  differen¬ 
tial  mobility  at  high  electric  fields  [2-3] .  In  addition, 
since  CdSnP2  has  a  bandgap  and  lattice  constant  close  to 
those  of  InP  it  seems  to  be  a  good  material  for  hetero¬ 
structure  devices  with  InP.  Growth  of  CdSnP2  as  platelets  [4] 
and  epitaxially  on  InP  [5]  has  so  far  yielded  material  with 
high  carrier  concentrations  and  low  mobilities  compared  to 
the  III-V  compounds.  In  fact,  attempts  to  produce  epitaxial 
layers  of  CdSnP2  on  InP  have  yielded  layers  with  very  poor 
morphology  and  much  lower  mobilities  compared  to  the 
platelets.  One  objective  of  this  project  was  to  grow  CdSnP2 
platelets  and  make  measurements  of  their  electrical  charac¬ 
teristics  and  to  improve  the  quality  of  epitaxial  layers 
grown  on  InP  in  an  attempt  to  grow  device  quality  material. 

As  in  the  CdSnP2/InP  system  ZnSnP2  growth  on  GaAs  is  of 
interest  for  the  study  of  heterostructure  devices.  In  addition 
ZnSnP2  is  an  attractive  material  for  the  study  of  chalcopyrite 
growth  on  the  III-V  sphalerite  structure  of  GaAs  for  numerous 
reasons  although  no  epitaxial  growth  experiments  have  yet 
been  reported.  The  lattice  constant  of  ZnSnP2  mismatches  that 
of  GaAs  by  only  0.04%  so  liquid  phase  growth  can  be  achieved 
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without  alloying  as  was  required  in  the  CdSnP2/InP  system.  Be¬ 
cause  the  c/a  ratio  of  the  tetragonal  ZnSnP2  structure  is  2.00 
growth  experiments  can  be  performed  on  various  substrate  orienta 
tions  without  incorporating  any  c-to-2a  mismatch.  This  allows 
us  to  study  the  effects  of  substrate  orientation  on  LPE  growth 
which  was  not  possible  in  the  CdSnP2/lnP  system.  The  second 
objective  of  this  project  was  to  study  the  LPE  growth  of  ZnSnP2 
on  GaAs  to  determine  if  ZnSnP2  is  a  potentially  useful  semi¬ 
conductor  material. 

CdSnP2  and  ZnSnP2  are  crystals  which  grow  with  the  chalco- 
pyrite  lattice  structure  shown  in  Figure  1.  This  structure  re¬ 
sembles  the  sphalerite  structure  of  the  III-V  compounds  with 
the  Column  III  element  sublattice  being  replaced  by  a  Column  II 
element  sublattice  and  a  Column  IV  element  sublattice.  CdSnP2 
is  a  direct  analog  of  InP  since  the  In  sublattice  is  replaced 
by  Cd  and  Sn  sublattices  where  Cd  and  Sn  straddle  In  on  the  peri 
odic  table.  To  accomodate  these  two  sublattices  the  size  of 
the  cubic  unit  cell  must  be  doubled.  Most  chalcopyrites  exhibit 
a  uniaxial  compression  along  their  longer  dimension  (ZnSnP2  is 
an  exception)  so  that  two  different  lattice  constants  are  needed 
to  describe  the  cubic  unit  cell:  a  and  c  as  shown  in  Figure  1 
where  c/a  £  2  for  all  chalcopyrites. 

Various  material  parameters  for  CdSnP2,  ZnSnP2,  GaAs, 
and  InP  are  presented  in  Table  1  [6] .  The  mismatch  of  CdSnP2 
on  InP  is  0.53%,  which  is  considered  small  for  some  growth 

systems  but  which  is  large  for  LPE  growth.  This  is  believed 
to  be  the  major  reason  for  poor  morphology  and 
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Figure  1 


Cubic  unit  cells  of  the  InP  sphalerite  structure 
and  the  CdSnP2  chalcopyrite  structure. 


poor  electrical  properties  of  CdSnP2  epitaxial  layers  on 
InP.  The  a  lattice  constant  of  ZnSnP2  is  a  good  deal  smal¬ 
ler  than  that  of  InP  so  a  chalcopyrite  alloy  of  CdSnP2  and 
ZnSnP2  might  be  produced  which  can  lattice  match  InP.  This 
was  the  impetus  behind  the  experiments  attempting  to  grow 
the  lattice  matched  chalcopyrite  alloy,  Zn  Cd,  SnP,  on 
InP.  Assuming  Vegard's  law  is  approximately  obeyed,  the 
alloy  composition  which  lattice  matches  {100}  InP  is 

Zn0.13Cd0.87SnP2* 

It  should  be  noted  that  CdSnP2  has  a  melting  point 
much  smaller  than  those  of  most  semiconductor  materials  and 
so  the  growth  of  CdSnP2  is  restricted  to  comparatively  low 
temperatures.  In  addition,  loss  of  volatile  constituents 
from  the  melt,  especially  phosphorus,  will  be  appreciable 
near  the  melting  point. 

Because  the  melting  point  of  ZnSnP2  is  closar  to  those 

of  the  III-V  compounds  the  growth  of  ZnSnP2  is  not  as 

restricted.  However,  ZnSnP2  makes  a  phase  transition  to  a 

° 

disordered  sphalerite  phase  at  720  C  [6] .  This  phase  has 
been  observed  in  solution  grown  platelets  in  which  cooling 

O 

rates  greater  than  5  C/hr  were  used  [9] .  Thus  LPE  growth 
of  ZnSnP2  on  GaAs  is  restricted  at  least  to  temperatures 

o  O 

below  720  C  and  cooling  rates  less  than  5  C/hr  in  order  to 
grow  the  chalcopyrite  phase. 
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2.1  CdSnP2  PLATELETS 

2.1.1  Growth  Procedure 

CdSnP2  was  initially  grown  as  platelets  from  a  Sn 
solution  to  make  material  property  measurements  and  for 
possible  use  as  a  source  material  in  the  LPE  system. 
Following  a  procedure  described  by  Buehler  et  al  [4], 
elemental  Cd,  Sn  and  red  P^  were  placed  in  a  vitreous 
carbon  boat  and  sealed  in  an  eight  inch  25  mm  O.D. 
quartz  ampoule  evacuated  to  a  residual  pressure  of 
2  x  10-6  torr.  All  elemental  materials  used  were  99.999% 
pure.  The  quartz  ampoule  was  nominally  1  mm  thick  to 
withstand  the  phosphorus  pressures  produced  on  heating. 
The  carbon  boat,  quartz  ampoule  and  quartz  sealing  plug 
had  been  cleaned  in  aqua  regia  for  24  hours  and  baked  out 
under  vacuum  prior  to  loading.  The  ampoule  was  heated 
over  a  period  of  six  hours,  allowed  to  homogenize  for  one 
hour,  and  then  slowly  cooled  to  230°C.  A  schematic  of  a 
typical  sealed  ampoule  is  shown  in  Figure  2.  Data  for 
these  growth  runs  are  given  in  Table  2 . 

The  resulting  ingots  containing  CdSnP2  in  Sn  were 
immersed  in  Hg  and  heated  to  about  60 °C  to  dissolve  the 
Sn.  The  CdSnP2  crystals  were  separated  from  the  bulk 
of  the  Hg-Sn  solution  by  passing  it  through  a  gauze  sieve 
and  the  remaining  solution  was  removed  in  33%  HNO^.  The 
resulting  CdSnP2  material  was  typically  in  the  form  of 


Figure  2  Sealed  ampoule  system  used  for  growing  CdSnP-  platelets. 
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10  mm  x  4  mm  x  3  mm  odd  shaped  single  crystals  mixed 
with  an  abundance  of  smaller  crystals  (2  mm  x  2  mm  x  2  mm) . 

2.1.2  Sample  Preparation 

Experimental  samples  were  prepared  from  the  larger 

platelets  by  chemical-mechanical  polishing  with  a  1%  Br 

in  methyl  alcohol  solution.  The  resulting  samples  were 

100  ym  to  400  ym  thick  with  parallel  flat  surfaces  of 
2  2 

area  0.25  cm  to  0.50  cm  . 

2.1.3  Contacts 

Ohmic  contacts  were  made  to  the  CdSnP2  samples 
using  10  mil  and  20  mil  diameter  Sn  spheres  alloyed  at 
300°C  for  two  minutes  in  a  high  purity  hydrogen  atmos¬ 
phere.  These  contacts  were  found  to  have  a  contact  ^ 

-4  2 

resistance  of  1.9  x  10  n  cm  at  300  K.  Although 
these  contacts  adhered  well  to  the  sample  surface  at 
300  K  they  tended  to  break  off  of  the  samples  when 
immersed  in  liquid  nitrogen.  No  contact  resistance 
measurements  were  made  at  77  K. 

Gold  plated  onto  the  surface  of  several  staples 
produced  softly  rectifying  contacts.  The  contacts  were 
alloyed  for  twenty  seconds  at  300 °C  in  a  hydrogen  at¬ 
mosphere  to  produce  ohmic  contacts.  No  contact  resis¬ 
tance  measurements  were  made  on  these  contacts. 

A1  contacts  evaporated  onto  the  surface  of  several 
samples  gave  rectifying  contacts  with  various  reverse 
breakdown  voltages  ranging  from  one  to  five  volts. 
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The  best  barriers  were  obtained  by  evaporating  A1  onto 
the  entire  sample,  masking  contact  regions  with  apiezon 
wax ,  and  removing  the  excess  A1  in  HF.  The  contacts 
produced  in  this  manner  had  higher  breakdown  voltages 
than  those  produced  by  A1  evaporation  through  a  glass 
or  metal  mask. 

2.1.4  Hall  and  Resistivity  Measurements 

The  majority  of  CdSnP2  samples  were  used  for  van 
der  Pauw  and  Hall  measurements  at  300  K,  77  K,  and  at 
elevated  temperatures  (300  K  to  625  K) .  These  samples 
were  prepared  as  described  above  with  four  10  mil  dia¬ 
meter  Sn  contacts  alloyed  to  the  sample  surfaces.  Those 
samples  which  were  measured  at  elevated  temperatures 
were  encapsulated  with  pyrolytic  Si02  grown  at  300 °C. 

The  encapsulant  layers  were  0.20  ym  to  0.40  ym  thick. 

The  samples  measured  at  300  K  were  n  type  with  a 

17  -3 

mean  carrier  concentration  of  1  x  10  cm  and  a  mean 

2 

mobility  of  725  cm  /V- sec.  The  better  samples  had 
carrier  concentrations  in  the  range  of  5  x  1016  to 
8  x  10^®cm”^  with  mobilities  of  1400  to  2000  cm2/V-sec. 
When  these  samples  were  measured  at  77  K,  no  substantial 
changes  in  the  carrier  concentrations  or  mobilities 
from  the  300  K  values  were  observed.  300  K  and  77  K 
measurements  were  made  using  a  2.15  kG  permanent 
magnet. 


# 


High  temperature  Hall  measurements  were  made  in 
the  apparatus  shown  schematically  in  Figure  3.  These 
measurements  were  made  in  hydrogen  purified  by  a 
palladium  diffusion  cell  using  a  flow  rate  of  15  ml/4nin. 
A  3.5  kG  magnetic  field  produced  by  a  four- inch  elec¬ 
tromagnet  was  used  for  these  measurements. 

High  temperature  measurements  of  the  encapsulated 
samples  were  limited  to  temperatures  below  700  K  due 
to  irreversible  changes  in  the  samples.  The  results 
of  measurements  on  a  typical  sample  are  plotted  in 
Figure  4.  In  these  experiments  the  temperature  was 
increased  from  room  temperature  and  measurements  were 
taken  every  25  K.  The  circles  in  Figure  4  are  the 
data  for  the  initial  set  of  measurements  and  the  squares 
are  from  a  subsequent  set  of  measurements  which  were 
made  after  the  sample  had  returned  to  room  temperature. 
These  data  indicate  a  steady  increase  in  the  carrier 
concentration  and  decrease  in  the  mobility  with  in¬ 
creasing  temperature  up  to  600  K.  Above  600  K  the 
changes  accelerate  and  are  irreversible  as  can  be  seen 
from  the  lower  mobilities  and  higher  carrier  concen¬ 
trations  at  low  temperatures  in  the  latter  measure¬ 
ments  . 

The  irreversible  changes  are  due  to  loss  of 
volatile  constituents  from  the  samples;  probably 
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phosphorus  due  to  its  high  vapor  pressure.  This  is 
supported  by  the  observation  of  numerous  cracks  in  the 
encapsulating  layers  of  the  samples  after  measurements 
were  made.  The  above  data  indicate  that  depletion  of 
volatile  components  from  the  samples  causes  an  increased 
n  type  carrier  concentration  and  a  decreased  mobility. 
Thus,  to  produce  epitaxial  layers  of  high  quality,  empha¬ 
sis  should  be  placed,  on  maintaining  a  high  concentration 
of  the  volatile  component  in  the  growth  melts. 


2.2  CdSnP2  GROWN  ON  InP 

2.2.1  Sealed  Tube  Epitaxial  Growth 

Epitaxial  growth  experiments  were  initially  per¬ 
formed  using  a  tipping  boat  arrangement  described  by 
Shay  et  al  [5]  which  is  shown  schematically  in  Figure  5. 
In  these  experiments  a  carbon  boat  was  used  to  contain 
the  substrate  and  growth  melt.  The  growth  melt  was 
homogenized  in  a  separate  step  by  loading  the  cleaned 
carbon  boat  with  99.999%  Cd,  Sn,  and  red  P^  and  sealing 
the  boat  in  a  cleaned  quartz  ampoule  evacuated  \.o  a  re¬ 
sidual  pressure  of  2  x  10~^  torr.  The  starting  materials 
were  measured  to  give  a  molar  composition  of  94%  Sn, 

2%  Cd,  and  4%  P  with  a  total  mass  of  10.6  grams.  This 
ampoule  was  then  heated  to  600°C  over  a  period  of  six 
hours,  allowed  to  homogenize  for  an  hour,  and  quenched. 
The  quartz  ampoule  was  broken  open,  a  substrate  was  in¬ 
serted  into  the  carbon  boat,  and  the  boat  was  sealed 
in  another  quartz  ampoule  after  evacuating  to  2  x  10  6 
torr  and  backfilling  with  argon  to  620  torr.  The 
{100}  Fe-doped  InP  substrate  had  been  prepared  by 
chemical-mechanical  polishing  one  face  with  4%  bromine- 
methyl  alcohol  followed  by  cleaning  in  boiling  tri¬ 
chloroethylene,  acetone,  and  methyl  alcohol,  and  a 
four  minute  etch  in  2%  bromine-methyl  alcohol.  The 
loaded  ampoule  was  held  in  a  vertical  furnace  as  shown 
in  Figure  5  and  was  rapidly  heated  to  525°C  and  allowed 
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Figure  5  Schematic  of  sealed  tube  tipping  LPE  system 
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to  equilibrate  for  fifteen  minutes.  The  furnace  was 
then  cooled  to  510°C  and  rotated  180°  so  that  the  melt 
solution  would  flow  over  the  substrate  and  out  of  the 
hole  in  the  substrate  holder  as  the  furnace  was  cooled 
10°C/hour.  When  the  furnace  had  cooled  to  360°Cf  it 
was  tipped  back  to  its  original  position  and  the  ampoule 
was  quenched.  The  separate  homogenizing  run,  argon 
backfill,  and  the  baffle  on  the  substrate  holder  were 
all  used  to  decrease  the  transport  of  volatile  con¬ 
stituents  from  the  melt  to  the  substrate  before  melt 
contact  was  made.  These  constituents  were  found  to 
cause  severe  pitting  of  the  substrate. 

The  epitaxial  growth  produced  in  this  manner  was 
very  rough  with  uneven  interfaces.  A  typical  cross 
section  of  a  sample  is  shown  in  Figure  6 .  The  epitaxial 
layers  were  generally  100  ym  thick.  No  Hall  measure¬ 
ments  were  made  on  these  layers. 

Because  of  the  poor  layer  quality,  little  time  was 
spent  working  with  this  system.  Instead,  an  open  tube 
H2  flow  system  was  designed  and  constructed  for  epitax¬ 
ial  growth.  Such  a  system  has  several  advantages  over 
the  closed  tube  system  in  that  the  apparatus  and  ma¬ 
terials  can  be  additionally  purified,  an  in  situ  etch 
can  be  used,  less  materials  are  used  in  each  growth 
run,  and  the  growth  runs  are  more  easy  to  control  over 
shorter  times. 
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2.2.2  Open  Tube  Flow  System 

The  open  tube  LPE  system  shown  schematically  in 
Figure  7  was  designed  and  constructed.  This  system 
consists  of  a  high  purity  graphite  sliding  boat  with 
stationary  melt  wells  and  a  slider  with  0.015"  deep 
substrate  wells.  This  boat  sat  in  a  40  mm  O.D.  quartz 
furnace  tube  through  which  purified  hydrogen  flowed. 
Between  growth  runs  the  system  was  baked  out  for  four 
hours  at  600°C  with  a  flow  rate  of  120  ml/min. 

The  {100}  InP  substrates  used  in  these  experiments 
were  chemically-mechanically  polished  on  one  side  and 
their  backs  were  lapped  using  5  urn  and  2  ym  grit  to 
a  thickness  of  400  ym.  These  substrates  were  then 
cleaned  and  etched  as  described  above.  The  final  chem¬ 
ical  etch  removed  approximately  30  ym  of  material  from 
the  surface.  The  {111}  and  {211}  InP  substrates  were 
prepared  by  lappir.g  both  faces  to  produce  substrates 
400  ym  thick  with  parallel  flat  faces  followed  by  clean 
ing  in  boiling  organic  solvents  and  a  chemical  etch  for 
two  minutes  in  1%  bromine-methyl  alcohol.  Substrates 
prepared  in  this  manner  would  fill  the  substrate  wells 
so  that  most  or  all  of  the  growth  solution  could  be 
removed  after  growth.  All  substrates  used  in  these 
experiments  were  Fe-doped  semi-insulating  InP. 

The  LPE  boat  was  loaded  by  placing  the  substrates 


in  the  substrate  wells ,  positioning  them  between  the 


Figure  7  LPE  sliding  boat  system 
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melt  wells,  and  then  loading  the  melt  wells  with  growth 
materials  or  etch  materials.  The  boat  was  then  loaded 
into  the  furnace  tube  and  rapidly  heated  to  the  homo¬ 
genization  temperature.  The  melt  was  allowed  to  homo¬ 
genize  for  one  to  two  and  a  half  hours  and  then  cooled 
to  the  growth  temperature:  typically  20°C  below  the 
homogenization  temperature.  At  this  point  the  sub¬ 
strate  was  etched  and  moved  to  contact  the  growth  so¬ 
lution.  The  furnace  was  then  cooled  at  a  constant 
rate:  generally  5°C/hr  or  10°C/hr.  A  twenty  second 
etch  of  the  substrate  in  pure  Sn  removed  the  top  5  um 
of  the  surface  and  gave  a  much  better  surface  for 
epitaxial  growth  than  an  unetched  substrate.  Such  an 
etch  was  used  in  all  of  the  growth  runs  except  the 
initial  few. 

In  the  first  forty -three  growth  runs  the  growth 
melt  was  made  up  of  99.999%  purity  ground  red  P4, 
pellet  Cd,  and  Sn  shot  with  molar  composition  15%  P, 

7%  Cd,  and  78%  Sn.  The  homogenization  temperature 
used  was  540 °C  and  the  initial  growth  temperature 
varied  between  490°C  to  525°C  ,  followed  by  a  10°C/hr 
cooling  rate  for  four  to  ten  hours.  This  procedure 
was  discontinued  because  of  the  inability  to  consis¬ 
tently  saturate  the  solution  at  a  given  temperature. 
This  was  due  to  loss  of  material  from  the  melts  before 
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the  melts  could  be  homogenized  and  during  growth. 
Epitaxial  material  was  grown  in  this  system  but  this 
material  was  in  the  form  of  isolated  single  crystals 
and  rough  epitaxial  layers.  The  best  CdSnP2  layer 
grown  in  this  system  is  shown  in  Figure  8. 

To  better  control  the  CdSnP2  concentration  in  the 
melts  so  that  the  saturation  temperature  could  be  more 
accurately  controlled,  two  changes  in  the  procedure 
were  made:  First,  a  Cd3P2  source  material  was  used  in 
place  of  Cd  and  P^  to  allow  for  more  complete  dissolu¬ 
tion  of  the  source  materials  in  the  Sn  solution.  Second, 
the  composition  of  Cd  and  P  in  the  melt  was  reduced  to 
lower  the  homogenization  and  initial  growth  temperatures 
and  reduce  the  loss  of  volatile  constituents  from  the 
melt. 

The  Cd3P2  source  material  was  grown  in  a  sealed 
tube  system  shown  schematically  in  Figure  9 .  In  this 
system  a  25  mm  O.D.  20  inch  quartz  ampoule  was  loaded 
with  a  boat  of  Cd  pellets  at  one  end,  red  P^  at  the 
other  end,  and  sealed  at  a  residual  pressure  of 
3  x  10~®  torr.  Stoichiometric  amounts  of  Cd  and  P, 
were  used  with  a  slight  amount  of  excess  P^  to  provide 
a  phosphorus  overpressure  in  the  ampoule.  A  total 
mass  of  30  grams  of  the  constituents  was  used  per 
growth  run.  This  ampoule  was  placed  in  a  furnace  and 


Figure  8  Cross-sectional  view  of  CdSnP2 
layer  grown  from  a  solution 
of  P*  and  Cd  in  Sn  on  {100}  InP 
after  chemical-mechanical  polish 
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heated  over  eight  hours  with  the  cadmium  at  750 °C  and 
the  phosphorus  at  a  mean  temperature  of  480 °C.  This 
cadmium  temperature  was  chosen  to  be  slightly  above  the 
melting  point  of  Cd3P2  and  the  phosphorus  temperature 
was  used  to  give  a  phosphorus  pressure  of  approximately 
two  atmospheres  in  the  ampoule.  The  ampoule  was  allowed 
to  homogenize  for  100  hours  and  was  then  cooled  to  room 
temperature.  Cd3P2  9rowt^  was  observed  near  the  center 
of  the  ampoule  in  a  region  where  the  temperature  had 
averaged  630 °C.  All  of  the  cadmium  and  phosphorus  had 
been  consumed.  Some  of  the  Cd.jP2  was  in  polycrystal¬ 
line  form  adhering  to  the  walls  of  the  ampoule  with 
large  single  crystals  growing  toward  the  center  of  the 
ampoule.  These  larger  crystals  were  used  for  the  source 
material. 

Liquid  phase  epitaxial  growth  was  resumed  using 
a  melt  solution  of  molar  composition  1.5%  Cd3P2  and 
98.5%  Sn.  The  homogenization  temperature  was  decreased 
to  the  430 °C  to  455 °C  range  with  a  two  and  a  half  hour 
homogenization  time  and  the  initial  growth  temperature 
was  typically  420°C.  Both  5°C/hour  and  10°C/hour  cooling 
rates  were  used  with  little  effect  on  the  growth. 

{100},  {111},  and  {211}  substrate  surfaces  were  used 
and  no  preference  for  any  orientation  was  observed. 

In  these  experiments  the  saturation  temperatures  were 
more  stable  although  there  were  still  some  fluctuations. 
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At  temperatures  above  455°C  excessive  Cd  and  P  loss  from 
the  melts  was  observed. 

The  layers  grown  in  this  manner  again  tended  to  have 
rough  surfaces  and  dendritic  growth.  A  typical  layer 
is  shown  in  Figure  10.  Hall  measurements  on  this  material 

indicated  that  it  was  n  type  with  a  carrier  concentration 
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of  1  x  10  cm  and  a  mobility  of  200  cm  /V-sec.  The 

high  temperature  Hall  measurements  indicate  that  the  poor 
material  quality  is  due  to  depletion  of  volatile  con¬ 
stituents  from  the  melt. 

X-ray  diffraction  measurements  were  made  on  some  of 
the  epitaxial  layers  and  some  of  the  data  are  plotted  in 
Figures  11  and  12.  The  measurements  indicated  a  mismatch 
of  -0.4%  to  -0.6%  as  compared  to  calculated  mismatches  of 
-1.9%  or  +0.5%  based  on  the  lattice  parameter  measurements 
of  Buehler  et  al  [4] .  The  differences  in  measured  and 
calculated  values  are  probably  due  to  differences  in  the 
growth  techniques.  The  peak  at  -1.8%  mismatch  in  Figure  12 
is  due  to  polycrystalline  regions  in  that  sample. 
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Figure  10.  Cross  section  view  of  CdSnP2 
layer  grown  from  a  solution 
of  Cd,P2  in  Sn  on  {211}  InP 
In  faeeT 


Figure 


2.3  Zn  Cd,  SnP,  GROWN  ON  InP 

We  believe  that  the  poor  quality  of  the  layers  grown 
with  the  previous  systems  can  be  attributed  primarily 
to  two  factors:  insufficient  lateral  mobility  of  the 
constituents  at  the  growth  surface  and  excessive  lattice 
mismatch.  The  lateral  mobility  of  the  constituents  can 
be  increased  by  increasing  their  thermal  energy:  that  is, 
by  increasing  the  initial  growth  temperature.  Because  of 
the  high  rate  of  evaporation  of  constituents  from  the 
melt,  several  alterations  were  made  to  the  system  for 
higher  temperature  work.  The  lattice  mismatch  can  be 
reduced  by  growing  a  suitable  alloy  of  CdSnP2  and  another 
chalcopyrite  so  that  the  alloy  lattice  constant  matches 
that  of  InP. 

To  minimize  the  number  of  system  parameters ,  chalco- 
pyrites  differing  from  CdSnP2  in  only  one  constituent 
were  reviewed  for  the  alloy  system.  ZnSnP2,  CdSiP2,  and 
CdGeP2  are  all  suitable  choices  for  this  purpose.  ZnSnP2 
was  chosen  since  the  composition  of  the  lattice  matched 
alloy  in  this  case  is  closer  to  CdSnP2  than  in  the  other 
alloys . 

Since  the  chalcopyrites  exhibit  uniaxial  compression 
it  is  not  possible  to  lattice  match  both  the  a  and  c  lattice 
constants  simultaneously.  Hence  all  growth  experiments 
were  performed  using  {100}  oriented  InP  substrates  so 
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that  the  a  lattice  constant  could  be  matched  with  the 
c  lattice  constant  extending  normal  to  the  growth  inter¬ 
face.  This  type  of  growth  could  be  anticipated  from  the 
x-ray  measurements  presented  previously.  Using  Vegard's 
law,  the  composition  of  the  alloy  which  should  lattice- 
match  InP  is  ZnQ  13CdQ  87SnP2. 

2.3.1  Revised  Open  Tube  System 

The  system  shown  schematically  in  Figure  13  was  used 
for  LPE  growth  of  ZnxCd1_xSnP2 .  This  apparatus  is  identi¬ 
cal  to  that  of  previous  systems  except  for  the  graphite 
boat  used.  In  the  new  boat  the  substrate  well  is  in  the 
base  with  the  melt  wells  extending  through  the  slider. 
During  the  melt  homogenization  process  the  slider  melt 
wells  are  aligned  with  those  in  the  top  section  and  the 
substrate  is  covered  by  the  graphite  slider.  This  reduces 
the  amount  of  volatile  constituents  transported  to  the 
substrat  and  reduces  substrate  pitting  due  to  reactions 
with  these  constituents.  To  maintain  a  saturated  solution 
during  homogenization,  excess  source  materials  were  added 
to  the  melt.  This  excess  material  floats  at  the  top  of 
the  solution.  To  initiate  growth  the  part  of  the  growth 
solution  contained  in  the  slider  is  separated  from  the 
rest  of  the  solution.  This  removes  the  excess  source 
material  before  the  melt  is  positioned  over  the  substrate. 
In  this  position  the  growth  solution  is  covered  by  the 
top  of  the  graphite  boat  to  reduce  evaporation  of  the 


Figure  13  Schematic  of  'JPE  boat  used  to  grow  Zn  Cd,  SnP 
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constituents.  By  using  the  apparatus  in  this  manner  the 
growth  solution  at  the  start  of  growth  is  saturated  and 
the  depletion  of  the  melts  by  evaporation  is  reduced. 

Because  the  vapor  pressure  of  phosphorus  is  much 
higher  than  that  of  the  other  melt  components,  it  is 
depleted  from  the  melt  faster.  To  compensate  for  this, 
it  was  decided  to  use  separate  sources  of  cadmium  and 
phosphorus  so  that  the  initial  phosphorus  concentration 
in  the  melt  could  be  increased.  For  this  reason,  elemental 
Cd  and  SnP3  grown  in  a  sealed  tube  system  to  be  described, 
were  chosen  as  sources  of  cadmium  and  phosphorus. 

2.3.2  Sealed  Tube  Synthesis  of  SnP3 

SnP3  was  grown  from  elemental  Sn  and  red  P4  using  the 

apparatus  shown  schematically  in  Figure  14.  It  consists  of 

a  20  mm  I. D.  x  25  mm  O.D.  x  9^  inch  ampoule  in  which  the 

Sn  and  P.  are  sealed  after  evacuation.  These  thick  wall 
4 

ampoules  were  required  to  withstand  the  high  phosphorus 
pressures  produced  during  growth.  Prior  to  loading,  the 
ampoules  and  sealing  plugs  were  cleaned  in  aqua  regia  and 
boiling  acetone  and  were  then  baked  out  under  vacuum. 
Stoichiometric  amounts  of  tin  and  phosphorus  were  used 
in  each  growth  run  with  a  few  grams  of  excess  phosphorus. 

A  total  mass  of  20  grams  were  used  per  growth  run.  The 
ampoules  were  loaded  and  sealed  at  a  residual  pressure 
of  4  x  10”6  torr.  The  ampoules  were  placed  into  the  fur¬ 
nace  so  that  the  Sn  was  in  a  6  inch  temperature  zone 
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flat  to  within  0.1°  C;  the  total  temperature  variation 
over  the  ampoule  being  less  than  2°  C.  The  Sn  was  spread 
over  a  length  of  approximately  4  inches  to  give  sufficient 
volume  for  the  growth  of  the  low  density  SnP^.  After 
loading  an  ampoule  into  the  furnace  it  was  heated  to 
500°  C  over  a  period  of  24  hours.  The  system  was  allowed 
to  homogenize  for  100  hours  and  then  was  cooled  to  room 
temperature . 

Upon  removal  from  the  furnace  the  ampoules  were  ob¬ 
served  to  contain  SnP,  and  excess  red  and  white  P..  The 

3  4 

ampoules  were  opened  in  an  argon  atmosphere  to  prevent 
combustion  of  the  phosphorus  and  the  SnP3  was  decanted. 
Prior  to  use  in  the  open  tube  epitaxial  system,  the  SnP^ 
was  ground  in  a  mortar  and  pestle  to  give  material  of  a 
more  convenient  size.  In  some  cases  the  SnP3  was  also 
cleaned  in  boiling  organic  solvents. 

The  stoichiometry  of  the  SnP3  was  calculated  from 
the  initial  mass  of  Sn  and  the  mass  of  the  product: 


m 


m 


M 


Sn 

SnP 

Sn 


=  Mass  of  Sn  loaded 

=  Mass  of  product, composition  SnPx 
x  A 

=  Atomic  mass  of  Sn  =  118-69  gms/mole 


Mp  =  Atomic  mass  of  P  =  30.9738  gms/mole 


x  =  3as  (  ^ .  ,) 

«P  \  "sn  / 
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The  values  of  x  calculated  in  this  manner  averaged 
around  2.8  and  were  slightly  smaller  than  the  theoretical 
value  due  to  loss  of  small  amounts  of  the  product  during 
decanting. 

2.3.3  Zn-Cd-SnP^-Sn  System 

Liquid  phase  epitaxial  growth  experiments  were  per¬ 
formed  to  grow  the  lattice-matched  chalcopyrite  alloy, 

Zn  Cd.  SnP-,  on  {100}  InP  using  the  previously  described 
apparatus.  Elemental  Zn,  Cd,  and  Sn,  all  of  5N  purity, 
and  the  synthesized  SnP3  were  used  as  source  materials  for 
the  system.  The  substrates  used  were  generally  Fe  doped 
semi-insulating  InP  prepared  as  discussed  previously  with 
a  few  experiments  performed  on  n  type  Sn  doped  substrates. 

A  typical  growth  run  consisted  of  a  one  hour  homo¬ 
genization  at  540°  C,  a  second  etch  of  the  substrate  in 
pure  Sn,  contact  of  the  growth  solution  with  the  substrate 
and  cooling  of  the  system  at  a  linear  rate  of  5°  C/hr. 

A  typical  melt  charge  consisted  of  0.2254  gm  of  Cd, 

0.3037  gm  of  SnP^,  1.0000  gm  of  Sn,  and  1.65  mg  of  Zn. 
The  equivalent  composition  of  such  a  melt  is  19.85% 

CdSnP2,  0.25%  ZnSnP2/  2.42%  excess  P,  and  77.48%  Sn. 

This  gives  a  ratio  of  ZnSnP2  to  CdSnP2  in  the  melt  of 
0.0125  which  is  about  one  tenth  of  the  calculated  ratio 
for  the  lattice  matched  solid.  All  growth  experiments 
were  carried  out  in  an  UHP  H2  atmosphere  with  a  flow  rate 
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through  the  tube  of  30  ml/min.  The  cooling  time  for 
these  experiments  varied  between  10  hours  and  22  hours. 

A  10  hour  cooling  time  gave  between  1  ym  and  2  ym  of 
growth  while  a  22  hour  tims  gave  approximately  15  ym  of 
growth  for  a  5°  C/hour  cooling  rate. 

A  typical  epitaxial  layer  is  shown  in  Figure  15. 

This  figure  shows  a  fairly  uniform  epitaxial  layer 
containing  numerous  pits  which  appear  as  dark  regions 
in  the  surface  view  and  as  notches  in  the  cross-sectional 
view.  These  pits  were  of  variable  depth  and  often 
extended  through  the  epitaxial  layer  to  the  interface. 

The  morphology  of  the  epitaxial  layers  is  very  de¬ 
pendent  upon  the  system  parameters.  For  an  initial  growth 
temperature  of  530°  C  the  growth  was  polycrystalline. 

A  10°  C/hour  cooling  rate  also  gave  polycrystalline 
material.  A  cooling  rate  of  4  0  C/hour  would  give  no 
growth  due  to  depletion  of  the  melt  by  evaporation.  The 
ZnSnP2  to  CdSnP2  ratio  in  the  melt  is  also  fairly 
critical  as  can  be  seen  in  Figure  16.  This  figure  shows 
cross  sections  of  the  epitaxial  layers:  one  grown  using 
the  parameters  given  above;  the  other  grown  from  a  melt 
with  a  6%  larger  Zn  mass.  The  latter  layer  is  poly¬ 
crystalline.  Values  of  growth  parameters  for  optimum 
epitaxial  growth  in  this  system  are  summarized  in 


Table  3. 


(b) 


Figure  15  (a)  surface  and  (b)  interface 

of  a  typical  Zn 
epitaxial  layer 
{100}  InP. 


xCdl  SnP2 

grown  on^ 
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Table  3  Typical  growth  parameters  for  the  Zn-Cd-SnP.-Sn 
system.  3 


SOURCE  MATERIALS:  Zn,  Cd,  SnP3,  Sn 

Melt  Composition 

CdSnP2  19.85  Molar  % 

ZnSnP2  0.25  Molar  % 

Sn  77.50  Molar  % 

P  2.40  Molar  % 

ZnSnP2/CdSnP2  =  0.0125 

Substrates:  {100}  Fe  Doped  InP 

Homogenization  Temperature:  540°C 

Cooling  Rate:  5° C/hour 

Growth  Duration:  10  to  22  hours 

in  aitu  Etch:  Pure  Sn 

Growth  Atmosphere:  UHP  H 2 

Flow  Rate:  30  ml/min 
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The  morphologies  of  epitaxial  layers  grown  in  this 
system  are  greatly  improved  over  those  of  the  previous 
systems  although  still  not  suitable  for  device  fabrication. 
The  high  segregation  coefficient  of  ZnSnP2  in  CdSnP2 
suggests  two  possible  impediments  to  growth  of  device 
quality  epitaxial  layers:  Zn  depletion  from  the  melt  due 
to  its  rapid  incorporation  into  the  growing  layer  and 
constitutional  supercooling.  Some  experiments  were 
performed  in  which  we  attempted  to  introduce  a  thermal 
gradient  in  the  melt  to  relieve  the  problem  of  constitu- 
tional  supercooling.  This  was  done  by  placing  the  boat 
on  a  sealed  chamber  through  which  N2  flowed  to  cool  the 
bottom  of  the  boat.  It  is  not  known  whether  or  not  this 
produced  any  thermal  gradient  in  the  melt.  However,  no 
substantial  improvement  in  the  epitaxial  layers  was  ob¬ 
served.  Another  impediment  to  smooth  layer  growth  may  be 
the  formation  of  antiphase  boundaries.  These  boundaries 
are  likely  to  form  since  the  epitaxial  material  has  a 
greater  number  of  sublattices  than  the  substrates. 

Possible  solutions  to  this  problem  may  be  to  grow  on 
substrates  oriented  slightly  off  of  the  <100>  direction 
or  to  grow  on  substrates  of  a  higher  index  orientation. 
These  experiments  have  not  been  performed. 


9 


-46- 


2.3.4  Hall  and  Resistivity  Measurements 

Van  der  Pauw  and  Hall  measurements  were  made  for  the 
layers  grown  on  semi-insulating  substrates  with  the  results 
shown  in  Figure  17.  All  epitaxial  layers  were  n  type. 

Most  had  room  temperature  carrier  concentrations  around 

18  “3  2 

10  cm  and  mobilities  around  100  cm  N  sec,  similar  to 

samples  grown  in  previous  systems.  The  resistivities 

_  t 

were  around  3  x  10  Q  cm.  There  are  no  substantial 

changes  in  these  values  for  measurements  made  at  77  K. 

2 

One  layer  had  a  mobility  of  2000  cm  /V-sec  with  a 

19-3  ... 

carrier  concentration  of  3  x  10  cm  .  The  resistivity 

-4 

of  this  sample  was  1  x  10  ft  cm.  This  sample  is  of 
interest  because  the  mobility  is  about  twice  that  which 
can  be  obtained  for  equivalent  III-V  compounds  at  similar 
carrier  concentrations.  This  sample  indicates  one  poten¬ 
tial  advantage  of  the  low  effective  mass  II-IV-V2  chalco- 
pyrites  over  the  III-V  sphalerite  semiconductors:  that  is, 
higher  mobilities  can  be  achieved  at  similar  carrier 
concentrations . 

The  solid  line  in  Figure  17  shows  the  theoretical 
dependence  of  the  mobility  on  the  carrier  concentration 
for  scattering  from  ionized  impurities  with  screening. 

This  curve  was  calculated  from  the  Brooks-Herring  model 
for  ionized  impurity  scattering  [7]  using  degenerate  sta¬ 
tistics  for  the  high  carrier  concentrations.  The  con¬ 
centration  of  ionized  impurities ,  N^,  in  this  model  was 


used  as  an  adjustable  parameter  which  was  set  so  that  the 

curve  matched  the  high  mobility  data  points.  The  value  of 

.  20  —3 

Nj  obtained  in  this  manner  was  Nj  =  1  x  10  cm  .  As  can 

be  seen,  there  is  reasonable  agreement  between  the  experi¬ 
mental  data  and  the  theoretical  curve.  This  indicates 
that  the  observed  decrease  in  mobility  with  decreasing 
carrier  concentration  is  mostly  due  to  close  compensation. 
Because  high  conductivity  inhomogeneities  in  an 

epitaxial  layer  can  give  anomalously  high  calculated 
mobilities  t 8 3  r  further  experiments  were  performed  on 
the  high  mobility  sample  to  determine  whether  the  mobility 
value  was  real.  First,  this  sample  was  cleaved  into  four 
pieces  and  the  pieces  were  remeasured.  The  effects  of  a 
conducting  inhomogeneity  on  the  resistivity  and  Hall 
measurements  can  usually  be  eliminated  with  this  procedure. 
The  mobilities  and  carrier  concentrations  of  each  piece, 
however,  were  very  close  to  those  of  the  original  sample. 
This  can  be  seen  from  the  tight  grouping  of  the  data  in 
the  high  concentration  region  of  Figure  17.  To  further 
substantiate  these  data,  measurements  of  the  Hall  constant 
as  a  function  of  magnetic  field  were  performed.  These 
results  are  plotted  in  Figure  18,  where  a  decrease  in  the 
Hall  constant  with  increasing  magnetic  field  is  shown. 

This  behavior  is  typical  of  homogeneous  n-type  material. 

If  there  had  been  a  conducting  inhomogeneity,  such  as 
those  which  give  an  anomalously  high  apparent  mobility, 


Figure  18  Measured  Hall  constant  as  a  function  of  magnetic  field. 
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the  Hall  constant  would  have  increased  with  increasing 
magnetic  field.  Thus,  the  mobility  and  carrier  concentra¬ 
tion  of  this  epitaxial  layer  are  apparently  not  due  to 
inhomogeneity  effects  and  appear  to  be  an  accurate  indica¬ 
tion  of  the  high  mobility  which  can  be  attained  in  this 
material. 

The  high  ionized  impurity  concentrations  in  these 
samples  is  probably  due  to  loss  of  volatile  constituents 
from  the  epitaxial  layers  during  growth.  This  is  supported 
by  the  results  of  the  high  temperature  Hall  measurements 
made  on  the  platelets,  which  showed  a  decreased  mobility 
and  increased  carrier  concentration  with  annealing  time. 
Thus,  annealing  the  samples  in  a  phosphorus  or  phosphorus 
and  cadmium  overpressure  should  give  increased  mobilities 
and  decreased  carrier  concentrations  as  is  reduced. 

Such  experiments  are  an  appropriate  subject  for  future 
investigation. 

2.4  ZnSnP2  GROWN  ON  GaAs 
2.4.1  Zn-SnP3~Sn  System 

Liquid  phase  epitaxial  growth  of  ZnSnP2  on  GaAs  was 
performed  in  the  revised  open  tube  system  described  in 
Section  4.1.  The  source  materials  used  were  5N  Sn  and  Zn 
and  SnP^  synthesized  as  described  in  Section  4.2.  Growth 
experiments  were  performed  using  initial  growth  temperatures 

O  O 

between  530  C  and  640  C  with  linear  cooling  rates  between 


1  C/hr  and  5  C/hr.  The  growth  periods  ranged  from  6  to 
18  hours. 

O 

For  cooling  rates  of  5  C/hr  the  grown  layers  were  thin 
(<5ym)  and  solid  material  was  observed  in  the  growh  solutions 
at  the  end  of  the  growth  runs.  This  suggests  that  this  cooling 
rate  is  too  large  to  allow  the  melt  constituents  to  propagate 
to  the  growing  surface,  resulting  in  supersaturation  and 

O 

crystal  nucleation  in  the  melt.  For  cooling  rates  of  1.0  C/hr 

O 

and  1.5  C/hr  the  layers  were  thicker  (5ym  to  35ym)  and  no 
solid  material  was  found  in  the  growth  solutions  at  the  ends 
of  the  runs . 

The  melts  used  for  the  growth  experiments  had  equivalent 
ZnSnP2  concentrations  ranging  from  1%  to  2%  with  slight  ($0.2%) 
excess  p  concentrations.  Figure  19  shows  the  dependence  of 
the  initial  growth  temperature  on  equivalent  ZnSnP2  concentration 
in  the  melt  for  saturated  or  slightly  supersaturated  melts. 

All  growth  runs  were  performed  in  an  UHP  H2  atmosphere  at  a 
flow  rate  of  30  ml/min.  At  this  flow  rate  and  temperature 
range  it  was  found  that  the  loss  of  phosphorus  from  the  melt 
was  not  nearly  as  much  of  a  problem  as  it  was  in  the 
ZnxCd^_xSnP2 (x=0 . 125)  growth  system.  Consequently,  the  growth 
experiments  were  much  more  reproducible  even  at  small  cooling 
rates.  Melt  homogenizations  of  1  to  2  hours  at  the  initial 
growth  temperature  were  used.  Between  growth  runs  the  system 

O 

was  baked  out  for  a  minimum  of  3  hours  at  850  C. 


EQUIVALENT  %  ZnSnP2  IN  MELT 


Figure  19  Initial  growth  temperature  as  a  function 
of  ZnSnP2  concentration  for  slightly 
supersaturated  solutions. 
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Initial  growth  experiments  were  performed  on  {100},  {110}, 
{111},  and  {211}  oriented  substrates  to  determine  the  optimum 
orientation  for  growth.  Typical  growth  on  these  orientations 
are  shown  in  Figures  20,  21,  and  22.  All  substrates  were  pre¬ 
pared  as  described  previously  except  for  the  final  chemical 
etch.  The  final  5:1:1  ’-^2^2  etch  cycle  consisted  of 

a  10  min.  cool-down  followed  by  a  2  to  4  min.  etch. 

{100},  {211},  and  {111}  Ga  oriented  substrates  produced 
growth  which  was  quite  rough  and  generally  unsuitable  for 
electrical  measurements.  {110}  and  {111}  As  oriented  substrates 
had  rough  interfaces  but  uniform  surfaces  with  very  good  mor¬ 
phologies.  This  is  consistent  with  the  observation  of  predomi¬ 
nant  {111}  facets  in  solution  grown  platelets  [6]  and  with 
experiments  in  which  we  have  observed  low  growth  rates  in  the 
{111}  direction.  Layers  grown  on  {110}  oriented  substrates 
revealed  much  better  interfaces  but  had  surfaces  with  more 
structure  and  tended  to  be  of  non-uniform  thickness.  Because 
of  this  the  majority  of  latter  growth  experiments  were  per¬ 
formed  using  {111}  As  oriented  substrates. 

All  layers  grown  on  {110}  and  {111}  As  oriented  substrates 
had  a  few  small  patches  of  Sn  adhering  to  their  surfaces.  These 
were  removed  by  ultrasonic  cleaning  in  a  solution  of  mercuric 
chloride  in  dimethyl  formamide.  This  is  in  contrast  to  the 
layers  of  Zn  Cd,  SnP-  grown  on  {100}  InP  and  the  layers  of 
ZnSnP2  grown  on  other  substrate  orientations  which  were  found 
to  be  completely  covered  with  a  Sn  layer  upon  removal  from 


the  LPE  boat. 


Figure  20  Surface  (220x)  and  cr< 
typical  ZnSnP-  growth 
(b)  {110}  oriented  Ga< 
is  4  ym  thick. 
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Figure  21  (a)  Surface  (220x)  and  cross-sectional  (450x)  views 

of  ZnSn?2  growth  on  {211}  Ga  oriented  GaAs. 

(b)  Surface  (llOx)  and  cross-sectional  (450x)  views 
of  a  40  yin  ZnSnP_  layer  grown  on  {211}  As  oriented  GaAs 


Figure  22  Surface  (220x)  and  cross-sectional  (450x)  views  of 

typical  ZnSnPj  layers  grown  on  (a)  {111}  Ga  oriented 
GaAs  and  (b)  fill}  As  oriented  GaAs.  Both  layers  are 
about  5  ym  thick.  The  structure  of  the  surface  view 
(b)  was  enhanced  using  a  Nomarski  interferometer. 
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For  {111}  As  oriented  substrates  the  best  growth  was 

O 

obtained  using  an  initial  growth  temperature  of  560  C  with  a 

O 

cooling  rate  of  1  C/hr,  the  smallest  cooling  rate  used  in  these 
experiments.  The  melts  for  these  growth  runs  had  an  equivalent 
ZnSn?2  concentration  of  1.4%  to  1.5%  with  one  gram  of  Sn 

O 

used  per  run.  For  initial  growth  temperatures  below  530  C 
no  layers  could  be  grown  and  for  initial  growth  temperatures 

O 

above  600  C  the  morphologies  and  electrical  characteristics  of 
the  layers  degraded  although  layers  could  be  grown.  This  may 
be  due  to  stoichiometric  defects  in  the  layers  probably  due  to 
phosphorus  loss.  The  grown  layers  were  generally  5  ym  to  15  ym 
thick  except  those  grown  at  high  temperatures  which  were  up 
to  30  ym  thick. 

2.4.2  Hall  and  Resistivity  Measurements 

Van  der  Pauw  resistivity  and  Hall  measurements  were  made  on 

the  grown  layers.  All  layers  were  found  to  be  p  type  with  carrier 

18  -3  21  -3 

concentrations  ranging  from  1x10  cm  to  1x10  cm  and  calcula- 

2  2 

ted  mobilities  ranging  from  10cm  /Vsec  to  70cm  /Vsec.  The  resis- 

4  —1 

tivities  of  these  layers  ranged  from  2.5x10  ftcm  to  1.2x10  ftcm. 

A  plot  of  carrier  concentration  as  a  function  of  inverse  initial 
growth  temperature  is  shown  in  Figure  23.  The  graph  reveals  a 
general  increase  in  carrier  concentration  with  increased  initial 
growth  temperature.  The  straight  line  interpolated  through  this 
data  suggests  that  the  stoichiometric  defects  which  are  thought 
to  cause  the  high  carrier  concentration  have  an  activation 
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IN VERSE  GROWTH  TEMPERATURE,  1/T  (K'1) 

Figure  23  Carrier  concentration  of 
epitaxial  layers  versus 
inverse  initial  growth 
temperature . 
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energy  of  about  5eV.  Only  slight  changes  in  these  data  were 
observed  between  300K  and  77K  measurements. 

A  few  experiments  were  performed  in  which  Column  VI  dopants 
were  added  to  the  growth  melts  in  an  attempt  to  make  the  layers 
n  type.  Addition  of  as  much  as  0.133  mg  of  Te  to  the  =  1  gm 
melts  had  no  noticeable  effect  on  the  electrical  characteris¬ 
tics.  This  may  be  due  to  a  large  segregation  coefficient  of 
Te  in  this  system.  Doping  experiments  are  an  appropriate  sub¬ 
ject  for  future  research  since  the  usefulness  of  ZnSnP2  as  an 
electronic  material  depends  on  the  ability  to  make  it  p  and 
n  type. 

2.4.3  X-Ray  Diffraction  Measurements 

X-ray  diffraction  measurements  were  made  on  the  layers 
grown  on  {111}  As  oriented  GaAs  using  a  double  crystal  diffrac¬ 
tometer  aligned  for  the  GaAs  333  reflection  for  Cu  Ka  radiation. 
Peak  widths  at  half  peak  intensity  of  10  seconds  of  arc  were 
reproducibly  obtained  with  {111}  GaAs  used  for  growth  sub¬ 
strates.  Some  results  of  these  measurements  are  shown  in 
Figures  24  and  25. 

Most  of  the  samples  measured  revealed  a  single  broadened 

GaAs  peak  with  some  slight  background  reflection  as  seen  in 

Figure  24.  This  suggests  that  the  epitaxial  material  lattice 

.  -3 

matches  GaAs  to  within  5x10  %  based  on  a  10  second  resolution. 

The  FWHM  peak  widths  ranged  from  45  seconds  to  85  seconds.  The 
low  intensity  background  reflection  may  be  due  to  the  rough 
interface  or  strain  incorporated  in  the  layers. 


Figure  24  Typical  double  crystal  x-ray  diffraction 
scan  for  ZnSnP2  on  {111}  GaAs. 
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Figure  25  Double  crystal  x-ray  diffraction 
scan  of  ZnSnP2  layer  grown  on 
{111}  GaAs  with  a  cooling  rate 
of  5° C/hr . 
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Some  of  the  thinner  layers  revealed  a  secondary  peak  at 

about  one  tenth  of  the  intensity  of  the  main  peak  separated 

by  4.9  minutes  of  arc  an  seen  in  Figure  25.  These  layers  were 

0 

3  ym  thick  and  were  grown  using  a  5  C/hr  cooling  rate.  The 
layers  had  much  smoother  interfaces  than  those  grown  using 
smaller  cooling  rates  and  no  broad  background  x-ray  intensity 
was  observed.  The  lattice  constant  of  the  layers  calculated 

O 

from  these  data  is  5.645  A  which  is  0.1%  smaller  than  the 
5.651  A  lattice  constant  observed  by  other  researchers  on 
solution  grown  platelets  [6,9].  The  FWHM  peak  widths  of  the 
main  peaks  were  20  to  30  seconds  and  those  of  the  secondary 
peaks  were  around  120  seconds. 

2.5  CONCLUSIONS 

Several  methods  for  the  growth  of  CdSnP2  platelets  and 
epitaxial  layers  on  InP  have  been  presented.  Work  with  these 
systems  suggests  that  the  lattice  mismatch  of  CdSnP2  on  InP 
is  a  significant  impediment  to  the  liquid  phase  growth  of 
device  quality  epitaxial  layers. 

Liquid  phase  growth  of  the  chalcopyrite  alloy, 

ZnCd1-xSnP2  was  investigated  as  a  method  for  obtaining  device 
quality  epitaxial  layers.  This  alloy  should  have  physical 
properties  close  to  those  of  CdSnP2  with  the  advantage  of 
being  lattice  matched  to  (100)  InP.  An  open-tube  LPE  system 
was  used  to  obtain  epitaxial  layers  of  the  alloy  which  were 
smoother  than  the  CdSnP2  grown  with  previous  systems.  The 
morphology  of  the  Zn  Cd,  SnP-  layers,  however,  was  still  not 
of  device  quality.  A  sample  having  a  mobility  of 


2  19  —3 

2,000  cm  /Vsec  with  a  carrier  concentration  of  3x10  cm 

was  produced  using  this  system.  This  suggests  that  it  may  be 

possible  to  effectively  utilize  the  low  electron  effective 

mass  of  this  material. 

A  method  for  the  growth  of  ZnSn?2  epitaxial  layers  on 
GaAs  has  been  discussed.  This  growth  was  studied  due  to  the 
small  lattice  mismatch  between  the  materials  and  because  the 
ZnSnP2  structure  is  amenable  to  growth  on  high  index  orienta¬ 
tions.  The  optimum  orientations  for  growth  were  found  to  be 
{110}  and  {111}  As.  The  ZnSnP2  layers  had  improved  morpholo¬ 
gies  over  the  CdSnP2  and  ZrfxCd^_xSnP2  (x^O.lZS)  growth. 


« 
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3.  VAPOR  PHASE  EPITAXIAL  GROWTH  OF  ZnGeAs2 

ZnGeAs2  crystallizes  in  the  chalcopyrite  structure  which 
has  tetragonal  symmetry.  It  can  be  viewed  as  two  stacked 
zincblende  unit  cells  as  illustrated  in  Figure  1.  ZnGeAs2 
may  be  considered  as  an  isoelectronic  analog  of  GaAs  in  which 
the  Ga  is  replaced  by  equal  proportions  of  Zn  and  Ge  to  give, 
on  an  average,  the  same  total  number  of  bonding  electrons  per 

O 

atom.  In  addition,  the  a  lattice  constant  of  ZnGeAs2  (5.671A) 
is  close  to  that  of  GaAs,  and  epitaxial  growth  of  ZnGeAs2  on 
{100}  GaAs  should  produce  a  lattice  mismatch  of  about  0.32%. 
The  close  structural  similarity  and  small  lattice  mismatch 
suggest  that  the  growth  of  epitaxial  layers  of  ZnGeAs2  on 
GaAs  surfaces  is  feasible. 

ZnGeAs2  has  previously  been  grown  in  bulk  form.  The 
large  solid  solubility  of  Ge  in  this  material  and  the  high 
mel'ing  point  of  Ge  tend  to  preclude  the  use  of  liquid  epitaxy 
to  grow  ZnGeAs2.  Vapor  phase  techniques,  however,  have  been 
applied  with  great  success  to  the  epitaxial  growth  of  III-V 
compounds.  The  objective  of  this  part  of  the  work  is  to 
attempt  to  grow  ZnGeAs2  epitaxially  on  GAs  with  a  VPE  open 
tube  flow  system. 


5.670  11.153  1.967 

5.672  11.153  1.966 


Figure  1  Crystal  structures  and  physical  properties  of 
GaAs  and  ZnGeAS2> 
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3.1  Zn-Ge-AsCl3  METHOD 

ZnGeAs2  is  one  of  the  least  studied  II-IV-V2  chalcopyrite 
compounds.  Vapor  phase  epitaxial  growth  of  ZnGeAs2  was  inves¬ 
tigated  for  the  first  time  with  a  Zn-Ge-AsCl^  system  [2] .  It 
was  an  open  flow  system.  The  reactor  is  shown  schematically 
in  Figure  2.  The  Zn  is  transported  by  passing  H2  over  it.  The 
AsC13  is  reduced  to  As4  and  HC1  upon  entering  the  hot  furnace, 
and  the  HC1  reacts  with  solid  Ge  to  produce  GeCl2.  At  lower 
temperatures  the  GeC^  will  undergo  a  disproportionation  reac¬ 
tion  to  transport  the  Ge.  Epitaxial  layers  with  varying  quanti¬ 
ties  of  Zn,  Ge,  and  As  were  obtained.  None  of  these  layers  had 
the  proper  stoichiometric  composition  and  were  apparently  Ge 
epitaxial  layers  with  high  doping  concentrations  of  Zn  and  As. 

Because  of  the  problems  associated  with  this  method,  other 
growth  processes  were  investigated. 

3.2  Zn-GeCl4-Ge-As4  METHOD 

Selection  of  a  chemical  process  for  epitaxial  growth  is  based 
on  a  number  of  factors,  including  availability  of  suitably  pure 
reactants,  their  compatibility  with  the  growth  apparatus,  and  the 
thermodynamics  of  the  process.  A  chemical  analysis  is  essential 
not  only  to  arrive  at  the  optimum  growth  conditions  but  also  for 
the  proper  apparatus  design.  In  order  to  qualitatively  predict  the 
direction  of  a  chemical  reaction  as  a  function  of  temperature,  it 
is  necessary  to  know  the  Gibbs  free  energy  change  AG  for  the  reaction. 
The  Gibb  free  energy  change  AG  for  a  reaction  is  related  to  the 
equilibrium  constant  K  by, 

AG  =  -  RT  In  K  =  -  4.574T  log  K  [3]. 


In  the  usual  formulation  of  a  chemical  reaction  if  AG  is 


negative,  the  reaction  proceeds  to  the  right  and  K  >  1 
while  if  AG  is  positive  the  reaction  proceeds 
to  the  left  and  K  <  1.  To  determine  a  suitable  reactor 
system,  the  chemical  equilibrium  constants  for  the  reac 
tions  of  Ge  in  a  chloride  system  [2,4]  are  shown  in 
Figure  3 . 

In  the  Zn-GeCl4-Ge-As4  system,  the  GeCl4  serves  as 
a  source  of  Ge  and  CI2  for  Ge  and  Zn  transport.  The 
GeCl^  is  kept  in  an  external  bubbler.  Since  GeCl4  is 
a  volatile  liquid  at  room  temperature,  palladium  puri¬ 
fied  hydrogen  is  bubbled  through  the  GeCl4  to  saturate 
it  with  GeCl4  vapor.  This  forms  the  input  to  the  part 
of  the  reactor  containing  the  Ge  boat.  The  primary 
reaction  that  takes  place  is, 

GeCl4  +  Ge  -**  2  GeCl2. 

The  equilibrium  constant  K  for  this  reaction  is  given 
[7]  by  the  equation. 


where  T  is  temperature  in  degrees  Kelvin,  and  the  vapor 
pressures  used  are  in  atmospheres.  The  GeCl2  serves  as 
a  means  of  transport  for  the  Ge . 


Hydrogen  is  passed  over  the  As^  boat  to  transport 
As4  vapor.  The  vapor  pressure  is  determined  by  the  tem¬ 
perature  of  the  As4  boat. 

The  Zn  vapor  reacts  with  the  GeCl2  vapor  to  deposit 
solid  Ge  and  form  ZnCl2  vapor  according  to  the  reaction, 

GeCl2  +  Zn  -*•  Ge  +  ZnCl2. 

The  equilibrium  constant  K  for  this  reaction  is  given  [6] 
by  the  equations. 


log  K  =  — 


4.574 


-52536.91 


-  1.98  log  T  +  42.91288 


and 


K  = 


P(ZnCl2) 


P(GeCl2)P(Zn)  ' 


where  T  is  in  degrees  Kelvin  and  the  vapor  pressures  are 
in  atmospheres. 

By  applying  this  chemical  analysis  and  mass  conser¬ 
vation  laws,  reactant  partial  pressures  and  molar  ratios 
for  a  given  set  of  growth  parameters  are  calculated.  It 
is  assumed  that  the  total  pressure  is  independent  of  tem¬ 
perature  and  constant  throughout  the  system  at  760  torr. 
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3.2.1  Growth  Experiments 

Attempts  to  synthesize  ZnGeAS2  were  carried  out  by 
employing  high  purity  Zn,  Ge,  As4,  and  GeCl4.  Figure  4 
shows  the  reactor  which  was  used  for  this  purpose.  It 
consists  of  a  quartz  tube  in  a  furnace  with  three  inde¬ 
pendently  controllable  zones.  The  left  hand  zone  controls 
the  As4  vapor  pressure,  the  center  zone  controls  the  Zn 
vapor  pressure,  and  the  right  hand  zone  is  the  reaction 
zone.  The  substrate  holder  and  boats  are  fabricated 
from  fused  quartz .  The  ultra  high  purity  H2  is  obtained 
using  a  palladium  diffusion  cell  and  is  used  as  a  car¬ 
rier  gas  in  the  reaction  tube  to  move  the  reactants  into 
the  reaction  zone.  Each  H2  flow  is  controlled  by  inde¬ 
pendent  needle  valves  and  pneumatic  flow  controllers  and 
monitored  using  rotameters  with  pyrex  or  316  stainless 
steel  floats.  The  Zn  and  As4  sources  are  contained  in 
quartz  boats  and  inserted  in  the  furnace  tube.  The 
GeCl4  is  contained  in  a  bubbler  and  its  temperature  could 
be  varied  from  -20°  to  30°C.  A  hydrogen  dilution  line 
which  bypasses  the  GeCl4  bubbler  is  also  provided.  It 
is  therefore  possible  to  vary  the  mole  fraction  of  GeCl4 
flowing  into  the  reaction  tube  by  either  changing  the 
vapor  pressure  at  the  GeCl4  source  while  keeping  the 
total  flow  constant,  or  by  adding  pure  hydrogen  through 
the  dilution  line  at  a  fixed  GeCl4  source  temperature. 


Epitaxial  layers  were  grown  on  chromium-doped  GaAs 

semi-insulating  substrates,  oriented  2°  off  the  {100} 

planes.  The  measured  resistivity  for  these  substrates 

7 

was  typically  10  ohm-cm.  Before  preparing  the  sub¬ 
strates,  the  wafers  were  cleaved  into  the  desired  size. 
These  were  then  degreased  by  boiling  successively  in 
trichloroethylene,  acetone,  and  methanol.  The  samples 
were  then  dried  with  bibulous  paper  and  nitrogen  gas. 
Prior  to  insertion  of  a  sample  into  the  reactor,  a  chemi¬ 
cal  etch  was  performed  to  give  a  clean  surface.  An 
etchant  consisting  of  5  H2S04:1  H2C>2:1  H20  was  mixed 
and  allowed  to  cool  for  4  minutes.  The  sample  was  then 
etched  for  6  minutes  and  rinsed  in  deionized  water  while 
being  careful  not  to  expose  it  to  air.  After  rinsing, 
the  sample  was  again  dried  with  bibulous  paper.  The 
sample  was  then  immediately  placed  in  the  sample  holder 
and  inserted  in  the  furnace. 

Following  the  insertion  of  the  seed  crystal  into 
the  reactor,  the  entire  system  was  flushed  with  hydrogen 
at  a  high  flow  rate  for  30  minutes.  The  reaction  tube 
was  then  heated  to  the  growth  temperature  profile  and 
the  hydrogen  flow  rates  for  both  barrels  were  set  at 
100  ml/min.  The  whole  system  was  allowed  to  stabilize 


for  one  hour,  then  the  hydrogen  in  the  first  barrel  was 
diverted  through  the  GeCl4  to  start  growth.  There  was 
no  vapor  etching  of  the  sample  before  growth.  Growth 
was  allowed  to  continue  for  2  hours.  To  terminate  growth, 
the  hydrogen  was  diverted  around  the  GeCl4  and  the  hydro¬ 
gen  flow  rates  were  reduced  to  30  ml/min  for  both  barrels. 
The  furnace  was  turned  off  immediately  and  allowed  to 
cool  by  opening  it. 

To  try  to  find  the  experimental  conditions  under 
which  ZnGeAs2  could  be  grown,  the  Zn  and  As4  temperatures 
were  set  at  600°  and  580°C,  respectively.  These  temper¬ 
atures  provided  a  large  enough  flux  of  these  species  to 
essentially  flood  the  reaction  zone.  With  the  Ge  tempera¬ 
ture  set  at  560°C  to  produce  the  disproportionation 
reaction  with  GeCl4,  the  temperatures  of  GeCl4  and  the 
reaction  zone  were  systematically  varied.  Over  most  of 
the  range  of  these  experiments,  little  or  no  growth 
occurred  on  the  {100}  GaAs  substrates.  With  substrate 
temperatures  around  720°C  and  GeCl4  temperatures  in  the 
range  from  4  to  10°C,  a  number  of  good  quality  epitaxial 
layers  were  obtained. 

3.2.2  Material  Properties 

The  properties  of  the  layers  grown  with  this  method 
were  established  by  analysis  of  surface  morphology  and 
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interfacial  characteristics.  Hall  and  resistivity  measure¬ 
ments,  and  x-ray  diffraction. 

Surface  morphologies  and  interfacial  characteristics  of 
epitaxial  layers  were  studied  with  a  metallurgical  micro¬ 
scope.  The  back  surface  of  the  sample  was  lapped  with 
5  and  2  vtm  grits  to  remove  the  epitaxial  growth  from  the 
back  of  the  sample  and  to  facilitate  cleaving  the  sample. 

The  interface  between  the  epitaxial  layer  and  substrate 
was  made  visible  by  a  delineation  etch.  The  delineation 
etch  consisted  of  1  HF:4  ^0:3  HNO^.  Figure  5  shows  the 
upper  surface  and  interface  of  one  of  the  epitaxial  layers. 
The  layers  have  a  very  high  degree  of  preferred  orienta¬ 
tion  and  their  interfaces  are  quite  flat  and  smooth. 

Ohmic  contacts  were  made  to  the  sample  using  20  mil 
diameter  Sn  spheres  alloyed  at  220°C  for  two  minutes  in 
a  high  purity  hydrogen  atmosphere.  Van  der  Pauw  Hall 
and  resistivity  measurements  were  then  used  to  determine 
the  resistivity,  carrier  concentration,  and  Hall  mobility 
of  the  grown  layers.  The  results  indicate  that  these 

epitaxial  layers  were  all  n-type.  Typical  values  of 

.  19  -3 

electron  concentration  and  mobility  were  5x10  cm  and 

2  + 

100  cm  /V-sec,  respectively.  Some  of  these  n  layers 


(a)  Surface 


I 


(b)  Interface 


Figure  5  Epitaxial  layer  grown  in  the 
Zn-GeCl^-Ge-As^  system. 
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were  also  grown  on  1015cm”^  n-type  GaAs.  The  specific 

contact  resistance  of  these  n-n+  heterostructures  was 
2 

about  0.1  ft  cm  . 

A  Geiger  diffractometer  was  used  to  examine  material 
deposited  in  regions  near  the  seedholder  in  the  reactor. 
This  material  was  prepared  for  diffraction  analysis  by 
grinding  it  into  a  fine  powder  and  placing  it  on  a  flat 
plate.  This  plate  was  then  inserted  into  the  specimen 
holder  of  the  diffractometer.  The  materials  were  identi¬ 
fied  from  their  diffraction  patterns  with  the  aid  of 
a  Hanawalt  index  and  file.  Analysis  of  this  material 
indicated  that  the  epitaxial  layers  grown  in  the  same 
region  were  mostly  germanium.  Figure  6  shows  a  diffrac¬ 
tion  pattern  of  these  Ge  deposits. 

3.2.3  Discussion 

Although  these  Ge  layers  are  quite  good,  we  were 
not  able  to  grow  the  compound  ZnGeAs2  under  these  con¬ 
ditions.  The  preparation  of  this  compound  appears  to  be 
very  critically  dependent  upon  the  reaction  conditions. 
Therefore,  several  growth  runs  were  made  with  different 
growth  parameters,  and  x-ray  diffraction  measurements 
were  made  to  determine  the  composition  of  the  material 
deposited  at  various  places  in  the  reactor.  Some  of 
these  results  are  shown  in  Table  2. 


Ge  deposits. 


Table  2  An  investigation  with  x-ray  diffraction  . 


600  5.3  500  29.84  0.48  35.47  0.0006  660 
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Typically,  the  factors  which  affected  the  nature 
of  the  deposits  most  were  substrate  temperature  and  reac¬ 
tant  ratio.  Polycrystalline  As4  was  formed  at  lower 
temperatures  while  Ge  and  Zn3As2  were  formed  at  high 
temperatures.  The  Ge/H2  ratio  or  Ge  input  flux  varia¬ 
tions  were  achieved  by  changing  the  vapor  pressure  at 
the  GeCl^  source  while  keeping  the  total  flow  constant. 
The  data  show  a  strong  Ge  growth  rate  dependence  with 
Ge/H2  ratio.  The  choice  of  reactant  vapor  pressure 
ratio  Zn:GeCl4:As4  appears  to  be  important.  When  the 
molar  ratio  of  GeCl2  in  the  reactant  mixture  is  high, 
GeCl2  is  further  transported  down  the  reactor  tube,  and 
Ge  is  deposited  out  by  reaction  of  GeCl2  with  H2 .  When 
the  molar  ratio  of  GeCl2  in  the  mixture  is  low,  large 
quantities  of  Zn3As2  platelets  form  inside  the  Zn  boat 
together  with  Ge  platelets. 

In  summary,  the  preparation  of  ZnGeAs2  with  the 
Zn-GeCl4~Ge-As4  system  was  found  to  be  much  too  strongly 
dependent  upon  the  various  growth  parameters.  When 
As4  was  transported  over  the  Zn  boat,  Zn3As2  platelets 
were  formed.  Apparently,  this  greatly  reduces  further 
transport  of  Zn  or  As4  vapor  and  results  in  the  depletion 
of  Zn  or  As4  in  the  reaction  zone.  This  appears  to  be 
the  main  problem  with  this  Zn-GeCl  -Ge-As.  system: 
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it  gives  excellent  control  of  the  Ge  flux,  but  very 
little  control  for  Zn  and  As^.  Previous  to  this  study, 
we  investigated  a  Zn-Ge-AsCl^  reactor  and  found  that 
it  produced  good  As4  control.  For  these  reasons,  we 
attempted  to  synthesize  the  compound  in  a  Zn-GeCl4-AsCl3 
reactor,  which  should  incorporate  the  best  features  of 
the  two  previous  systems  without  their  disadvantages . 

3.3  Zn-GeCl4-AsCl3  METHOD 

In  the  Zn-GeCl4~AsCl3  system,  the  AsC13  serves  as  a 
source  of  As4  vapor  while  the  GeCl4  serves  as  a  source  of 
Ge.  The  AsC13  is  kept  in  an  external  bubbler.  Palladium 
purified  hydrogen  is  bubbled  through  the  AsC13  and  the 
AsC13  saturated  vapor  is  passed  into  the  first  barrel  of 
the  reaction  tube.  The  initial  reaction  taking  place  when 
the  AsCl3  and  H2  gas  mixture  enters  the  hot  reaction  tube  is, 

4  AsC13  +  6  H2  ->  As4  +  12  HC1. 

Since  the  equilibrium  constant  is  very  large,  this  reaction 
is  assumed  to  go  to  completion.  Therefore,  the  input 
pressures  for  As4  and  HC1  are  computed  as, 

P(As4)  =  1/4  P(AsC13), 

and 

P(HC1)  =  3  P  (AsC13) . 

These  input  pressures  are  then  normalized  for  a  total 
system  pressure  of  760  torr. 


In  the  second  barrel  of  the  reactor  system,  hydro¬ 
gen  is  passed  over  a  Zn  boat  to  transport  elemental  Zn 
vapor.  The  vapor  pressure  of  the  Zn  is  determined  by 
the  temperature  of  the  Zn  boat. 

The  GeCl4  is  kept  in  an  external  bubbler.  Hydrogen 
is  bubbled  through  the  GeCl^  to  saturate  it  with  GeCl4 
vapor.  Various  reactions  of  Ge  in  a  chloride  system 
were  examined  carefully  (Figure  3).  The  two  chemical 
reactions  involved  in  the  third  barrel  of  this  reactor 
system  are, 

% 

GeCl4  +  2  H2  -*•  Ge  +  4  HC1, 

and 

Ge  +  2  HC1  -*■  GeCl2  +  H2- 

When  the  GeCl4  and  H2  gas  mixture  enters  the  hot  reactor 
tube,  the  GeCl4  readily  reacts  with  H2  at  the  temperature 
of  the  first  zone  and  deposits  out  Ge.  HC1  vapor  then 
reacts  with  solid  Ge  to  produce  GeCl2  vapor.  The  GeCl2 
vapor  then  serves  as  a  means  of  transport  for  the  Ge. 
Calculations  of  the  reactant  partial  pressure  and  mole 
ratio  are  based  on  the  above  analysis. 

3.3.1  Growth  Experiments 

The  growth  system  employs  Zn,  GeCl4  and  AsCl-j  as 
reactants  with  ultra  high  purity  H2  as  the  carrier  gas. 
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High  purity  elemental  Zn  (99.999%)  and  electronic  grade 
GeCl4  and  AsCl^  are  used  as  starting  materials.  The 
reactor  consists  of  a  44  mm  OD  quartz  tube  having  sepa¬ 
rate  reactant  inlets  and  a  furnace  with  three  indepen¬ 
dently  controllable  zones  (Figure  7) .  The  left  hand 
zone  controls  the  Zn  vapor  pressure,  while  the  right 
hand  zone  is  the  reaction  zone .  The  Zn  boat  and  sub¬ 
strate  holder  are  fabricated  from  fused  quartz.  The 
GeCl ,  and  AsCl,  enter  the  reactor  downstream  from  the 
Zn  source.  Both  GeCl^  and  AsCl^  are  contained  in  gas 
saturators  held  at  the  same  temperature  by  using  a  re¬ 
frigerated  circulator.  An  additional  H2  flow  is  also 
provided  to  dilute  the  GeCl4  if  necessary.  H2  flow  rates 
through  different  inlets  of  the  reactor  are  adjusted  by 
needle  valves  and  monitored  with  rotameters  specifically 
calibrated  for  H2» 

Attempts  were  made  to  grow  epitaxial  layers  on 
Cr-doped  semi- insulating  GaAs  substrates.  The  substrate 
is  etched  in  a  5  H2S04:1  H202:1  H20  solution  as  previ¬ 
ously  discussed.  After  etching,  the  substrate  is  placed 
on  the  sample  holder  and  inserted  in  the  furnace.  The 
reaction  tube  is  then  flushed  with  H2  and  heated  to  the 
growth  temperature  profile.  There  is  no  vapor  etching 
of  the  sample  before  growth.  The  reactant  flows  of 


Figure  7  The  Zn-GeCl  .-AsCl..  reactor  system. 
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GeCl4  and  AsCl^  are  initiated  and  terminated  simultane¬ 
ously.  Following  deposition  at  the  growth  temperature, 
the  reaction  tube  is  then  cooled  to  room  temperature 
rapidly  to  terminate  the  transport  of  Zn  and  minimize 
surface  erosion  of  the  sample.  A  systematic  study  of 
the  many  experimental  parameters  was  carried  out.  The 
growth  parameters  used  are  shown  in  Tables  3  and  4 ,  in 
which  we  specify  the  temperature  profile  in  which  growth 
was  achieved,  the  reactant  flows,  total  flow  rates, 
growth  composition,  and  other  parameters. 

3.3.2  Material  Properties 

The  properties  of  the  material  grown  by  this  process 
were  established  by  analysis  of  surface  morphology  and 
interfacial  characteristics,  Hall  and  resistivity  mea¬ 
surements,  Auger  electron  spectroscopy,  and  x-ray  dif¬ 
fraction.  This  analysis,  however,  indicated  that  these 
layers  were  either  Ge  or  Zn3As2. 

Figure  8  shows  the  upper  surface  and  interface  of 

one  of  the  Ge  layers  which  appears  to  be  polycrystalline. 

* 

A  significant  amount  of  Zn  diffusion  into  the  GaAs  sub¬ 
strate  is  revealed  in  the  cross-sectional  view.  Zn3As2 
layers  were  also  synthesized  in  this  system.  Figure  9 
shows  the  surface  morphology  of  a  Zn3As2  layer. 


Growth  parameters  for  the  Zn-Ge”l4-AsCl,  system 
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Table  4  Growth  parameters  for  the  Zn-GeCl^-AsClo  system 


(b)  interface 

Figure  8  Ge  layer  grown  in  the 
Zn-GeCl4-AsCl2  system. 


9  Surface  morphology  of  a 
Zn.jAS2  layer. 
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Localized  stresses  in  the  layer  results  from  large  lattice 
mismatch  between  the  Zn3As2  layer  and  GaAs  substrate. 

The  Ge  layers  on  Cr  doped  GaAs  substrates  were  typi¬ 
cally  p-type  with  carrier  concentrations  of  about 

19-3  2 

7x10  cm  and  hole  mobilities  of  40  cm  /V-sec.  The 

Zn3As2  layers  tend  to  crack  easily.  Therefore,  no  Hall 

and  resistivity  measurements  were  made  on  these  layers. 

The  composition  of  the  material  synthesized  was  con¬ 
firmed  by  Auger  electron  spectroscopy  and  x-ray  powder 
diffraction  analysis. 

Auger  electron  spectroscopy  was  performed  to  deter¬ 
mine  the  composition  of  the  sample  surfaces.  Figure  10 
shows  the  polycrystalline  surface  morphology  of  Sample 
J14.6.  The  sample  appears  to  have  a  solid,  although 
multi-crystalline  surface.  Auger  spectra  for  this  layer 
before  and  after  sputter  etching  are  shown  in  Figure  11. 

The  spectra  indicate  that  the  large  crystals  on  the  sur¬ 
face  are  germanium.  The  surface  morphology  of  Sample  J14.7 
(Figure  12)  shows  dendritic  structure  with  white  patches 
present  on  the  surface.  The  Auger  spectra  of  Sample  J14.7 
in  the  area  of  the  white  precipitate  are  shown  in  Figure  13. 


Figure  10  S-ylid  polycrystalline 

structure  of  Sample  J14.6. 
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Figure  11  Auger  spectra  for  the 

surface  of  Sample  J14.6, 
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Figure  13  Auger  spectra  for  Sample 
J14.7  in  the  area  of  the 
white  precipitate. 
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The  spectra  show  that  the  patches  of  lighter  material 
are  Znx,  0x,  and  Clx-  Both  samples  are  contaminated  by 
carbon  and  oxygen. 

The  Debye- Scherrer  powder  x-ray  diffraction  method 

was  employed  to  analyze  the  composition  of  the  material 

grown.  Cu  K  radiation  with  a  Ni  filter  was  used  as  the 

o 

x-ray  source  at  a  wavelength  of  1.5405  A.  Platelets  de¬ 
posited  in  the  region  near  the  substrate  in  the  reactor 
were  prepared  for  diffraction  analysis  by  grinding  them 
into  a  fine  powder.  The  powder  was  mixed  with  a  small 
amount  of  glue  and  then  rolled  or  extruded  into  a  small 
cylinder.  The  powder  specimen  was  then  inserted  into  the 
specimen  holder  of  the  Debye  camera.  Compounds  are  iden¬ 
tified  by  comparison  of  their  powder  x-ray  diffraction 
spectra  with  data  in  the  ASTM  x-ray  powder  data  file  [5] 
Several  powder  diffraction  patterns  were  obtained. 

The  powder  photograph  yields  a  set  of  numerical  values 
for  various  interplanar  spacing  of  the  crystal  sample. 

The  interplanar  spacing  of  the  sample  can  be  calculated 
from  the  Bragg  equation, 

d  "  2  "sin  $ 


where  d  is  the  interplanar  spacing,  X  is  the  wavelength 
of  the  x-ray  source,  and  8  is  the  angle. 
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Figure  14  shows  powder  diffraction  photographs  of  Samples 
J14.6  and  J23.5.  Table  5  shows  the  ASTM  powder  data  and 
d  values  of  Sample  J14.6  calculated  from  its  diffraction 
pattern.  Table  6  shows  the  ASTM  powder  data  for  Zn3As2 
and  d  values  of  Sample  J23.5  calculated  from  its  diffrac¬ 
tion  photograph.  Results  show  that  Sample  J14.6  is  Ge 
and  Sample  J23.5  is  Zn3As2. 

3.3.3  Discussion 

Ge  and  Zn^ASj  layers  have  been  synthesized  in  the 
Zn-GeCl4~AsCl3  system  with  different  growth  parameters. 
The  existence  of  these  materials  was  confirmed  with  Auger 
electron  spectroscopy  and  x-ray  powder  diffraction 
analysis.  A  transition  from  Ge  to  Zn3As2  with  growth 
conditions  was  observed  as  indicated  in  Table  4.  All 
tnese  layers  were  synthesized  with  the  same  growth  tem¬ 
perature  profile  but  different  input  reactant  flow  rates. 
Also,  the  substrates  were  located  very  close  to  the 
three  barrel  section  of  the  reactor. 

The  growth  parameters  were  examined  to  see  what 
determines  whether  Ge  or  Zn3As2  will  be  deposited  on  the 
substrate.  Reactant  partial  pressures,  reactant  mole 
ratios,  and  the  growth  composition  of  the  samples  are 
summarized  in  Table  7.  From  the  data  in  Table  7  it  was 
found  that;  (1)  considerable  excesses  of  Zn  and  HC1 
vapors  exist  in  the  growth  process;  (2)  Ge  will  deposit 
on  the  substrate  if  the  molar  content  of  GeCl2  in  the 


Table  5  Interplanar  spacing  values  for  Sample  J14.6 


(a)  ASTM  power  data  for  Ge 


(b)  Interplanar  spacing  values  for 
Sample  J14.6  calculated  from 
its  diffraction  photograph 
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Table  6  Interplanar  spacing  values  for  Sample  J23.5 


(a)  ASTM  powder  data  for  Zn2As2 


(b)  Interplanar  spacing  values 
for  Sample  J23.5  calculated 
from  its  diffraction  photo¬ 
graph 
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Table  7  Reactant  partial  pressures,  reactant  mole  ratios,  and  growth  composition 


J19.1  59.723  10.713  0.6496  30.761  658.15  1:0.179:0.010 
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reactant  mixture  is  high  enough.  In  addition,  large 

numbers  of  Ge  platelets  were  always  formed  at  the  outlet 
of  the  third  barrel  in  the  reactor  during  the  growth 

experiments.  The  equilibrium  constant  of  the  reaction, 

GeCl2  +  Zn  -*•  Ge  +  ZnCl2, 
was  then  examined. 

The  chemical  equilibrium  calculation  and  experimen¬ 
tal  results  indicated  that  GeCl2  vapor  readily  reacts 
with  Zn  vapor  to  deposit  solid  Ge  and  form  ZnCl2  in  the 
reaction  zone.  If  the  vapor  pressure  of  GeCl2  is  large 
enough,  Ge  will  also  deposit  on  the  substrate  which  is 
located  very  close  to  the  three  barrel  section  of  the 
reactor. 

The  premature  reaction  of  GeCl2  with  Zn  results  in 
the  depletion  of  GeCl2  vapor  in  the  reactant  gas  mixture 
and  this  is  the  explanation  for'  the  failure  to  synthe¬ 
size  ZnGeAs2  in  this  system. 
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3.4  PROPOSED  HYDRIDE  METHODS 

Two  serious  problems  are  found  in  the  chloride  methods 
of  synthesizing  ZnGeAs2:  (1)  premature  reactions  prevent 
further  transport  of  source  materials?  (2)  the  complexity 
of  the  chemical  reactions  increase  the  difficulties  in 
controlling  growth  parameters  and  synthesizing  the  desired 
compound.  However,  two  proposed  hydride  growth  systf- 
can  be  used  to  overcome  these  disadvantages. 

3.4.1  Zn-GeH4-AsH3  Method 

The  first  proposed  hydride  system  employs  elemental 
Zn,  GeH4  and  AsH3  as  reactants  with  H2  as  the  carrier  gas. 
The  reactor  system  is  shown  in  Figure  15.  The  reactor 
consists  of  a  quartz  tube  having  separate  reactant  inlets. 
The  Zn  is  transported  into  the  deposition  region  by  passing 
H2  gas  over  a  Zn  melt.  The  GeH4  and  AsH3  reactants  enter 
the  reactor  downstream  from  the  Zn  source,  and  the  overall 
growth  reaction  is  as  follows: 

Zn  +  GeH4  +  2  AsH3  -*•  ZnGeAs2  +  5H2. 

The  thermal  decomposition  of  GeH4  and  AsH3  was  chosen  for 
several  reasons: 

(1)  Both  GeH4  and  AsH3  are  commercially  available  in 
high  purity  and  are  stable  gases  at  room  temperature. 

(2)  Deposition  can  be  achieved  at  relatively  low 


temperatures . 
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(3)  There  are  no  reactive  by-products  of  the 
decomposition  reactions. 

3.4.2  Zn(CH3) 2-GeH4-AsH3  METHOD 

The  second  growth  system  proposed  to  synthesize 
ZnGeAs2  is  a  horizontal,  RF  induction  heated,  open  tube 
chemical  vapor  deposition  system.  The  reactants  are 
Zn(CH3)2,  GeH4,  and  AsH3  with  H2  as  the  carrier  gas. 

The  liquid  dimethylzinc  (DMZn)  is  kept  in  a  bubbler  and 
the  hydrogen  carrier  gas  is  bubbled  through  the  liquid 
to  saturate  and  transport  its  vapor.  The  deposition 
of  ZnGeAs2  involves  the  pyrolysis  of  these  reactants 
according  to  the  overall  reaction, 

Zn(CH3)2  +  GeH4  +  2  AsH3  -►  ZnGeAs2  +  2  CH4  +  4  H2 . 

The  advantage  of  using  these  reactants  is  that  all  three 
can  be  transported  in  the  gas  phase  to  the  deposition 
region.  Furthermore,  transport  of  Zn(CH3)2  vapor  is 
readily  and  almost  as  easily  controlled  as  the  gaseous 
reactants,  GeH4  and  AsH.^. 


3.5  CONCLUSIONS 


Two  chloride  systems,  Zn-GeCl4-Ge-As4 ,  and  Zn-GeCl4~ 
AsCl3  were  employed  for  the  epitaxial  growth  of  ZnGeAs2 
on  GaAs .  Since  the  preparation  of  ZnGeAs2  was  so  cri¬ 
tically  dependent  upon  various  growth  parameters,  there 
was  little  success.  With  careful  choice  of  reactant, 
the  Zn-GeH4-AsH3  and  Zn  (CH.^)  2-GeH4~AsH3  systems  could  be 
used  to  overcome  the  difficulties  encountered  in  the 
chloride  systems.  Both  systems  should  provide  easy 
transport  of  reactants  and  precise  control  of  reactant 
ratio.  It  is  believed  that  epitaxial  ZnGeAs2  could  be 
successfully  deposited  on  GaAs  with  either  of  these 
hydride  systems. 
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4 .  HETERO JUNCTION  THERMODYNAMICS 

The  recent  interest  in  heterostructures  for  lasers, 
transistor  emitters,  real-space  transfer  devices,  quantum 
well  lasers,  modulation-doped  superlattices,  and  other  de¬ 
vices,  has  prompted  a  re-examination  of  basic  hetero junction 
theory.  To  understand  quantitatively  the  operation  of 
these  devices  it  is  necessary  to  determine  how  the  discon¬ 
tinuity  in  the  energy  gaps  of  a  heterostructure  is  divided 
between  the  conduction  and  valence  bands: 

Aeg  -  Aec  +  Aev  (1) 

The  basic  assumption  of  heterostructure  theory  is  the  "elec¬ 
tron  affinity  rule":  that  is,  the  conduction  band  discon¬ 
tinuity  is  simply  the  difference  between  the  electron 
affinity  of  the  two  materials, 

Aec  =  xl  "  x2  '  (2) 

presumably  referenced  to  the  infinite  vacuum  level. 

4 . 1  CURRENT  PROBLEMS 

Kroemer  [1]  has  pointed  out  a  number  of  problems  with 
this  concept,  the  most  important  of  which  are:  (1)  assum¬ 
ing  the  rule  is  valid,  the  band  discontinuities  need  to  be 
known  to  an  accuracy  of  about  a  kT,  while  the  affinities, 
usually  referenced  to  a  near  vacuum  level,  cannot  be  meas¬ 
ured  that  accurately;  (2)  the  basic  differences  between 
free  surfaces  and  solid  interfaces  are  so  great  that  there 


is  no  fundamental  reason  for  expecting  any  relationship 
between  electron  affinities  and  band  discontinuities.  Thus, 
to  resolve  this  problem,  it  is  necessary  to  forget  about 
phenomenological  surface  properties  and  look  at  the  inter¬ 
face  problem  from  a  microscopic  point-of-view. 

In  spite  of  these  arguments,  Shay  et  al  [2]  compared 
measured  band  discontinuities  for  several  hetero junctions 
with  measured  electron  affinities  from  photoemission  data 
and  calculated  affinities  from  the  dielectric  two-band  model 
They  found  that  the  agreement  between  the  measured  discon¬ 
tinuities  and  those  calculated  from  the  affinity  rule  was 
very  good.  In  addition,  Frensley  and  Kroemer  [3 J  found  from 
a  pseudopotential  calculation  of  a  hetero junction  interface 
that  the  affinity  rule  is  approximately  correct.  Although 
this  is  a  somewhat  surprising  result,  Kroemer' s  original 
objections  are,  of  course,  still  valid. 

More  recently,  Adams  and  Nussbaum  1 4 J  have  examined 
the  affinity  rule  from  an  electrostatic  point-of-view.  They 
suggest  that  instead  of  keeping  the  vacuum  level  continuous 
(the  electron  affinity  rule)  it  may  be  more  appropriate  to 
assume  continuity  of  the  intrinsic  level.  This  paper  is 
interesting  because  it  also,  although  unintentionally,  dem¬ 
onstrates  another  problem  with  heterojunction  theory:  the 
usual  equilibrium  and  non-equilibrium  formulations  in  the 
depletion  approximation  are  not  self-consistent.  This 
inconsistency  occurs  regardless  of  which  basic  assumption 


is  used 
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The  problem  arises  in  the  following  manner:  For  sim¬ 
plicity  assume  an  abrupt  p-n  hetero junction.  In  equilibrium, 
the  built-in  potential  of  the  space-charge  region  is  given 
by. 


„  ,  kT  ,  (  Va  \ 

V0  ‘  T  ln  l  n~n7~  )  ' 
^  \  m  ip  / 


(3) 


where  Nd  and  N&  are  the  donor  and  acceptor  concentrations  on 

the  n  and  p  side  of  the  hetero junction,  respectively,  and 

n.  and  n.  are  the  intrinsic  concentrations.  When  the 
in  ip 

hetero junction  is  forward-biased  with  a  voltage  V,  in  the 
depletion  approximation  the  holes  injected  into  the  n-side 
at  the  edge  of  the  depletion  region  is  usually  found  to  be. 


P„  =  P„  exp 
n  p 


[^]  ■ 


(4) 


If  V  is  turned  off,  the  excess  holes  disappear  and 


pn  = 


n2 

in 

N. 


and 


V  =  —  In 
0  q 


P  =  N  , 
P  a 


/NdNa\ 

ter 


(5) 


Comparing  Eqs.  (3)  and  (5),  we  see  "hat  this  analysis  is 
only  valid  for  a  p~n  homo junction. 

The  problem  we  wish  to  consider  is  the  latter  one, 
although  the  electron  affinity  rule  is  a  natural  consequence 
of  our  analysis.  We  will  show  that  the  inconsistency  arises 
from  the  neglect  of  a  thermodynamic  force  term,  which  is  well 
understood  but,  apparently,  not  widely  known. 


4.2  THERMODYNAMIC  ANALYSIS 


For  any  solid  it  can  be  shown  [5]  that  the  entropy 
source  strength  (entropy  production  rate  per  unit  volume)  is, 


=  J. 


A. 

_1 

T 


v .  > 

J  ~ 


(6) 


In  this  equation  Jq  is  the  heat  flux  (flux  of  internal  ener- 
gv) ,  are  the  particle  fluxes  for  each  i  component,  are 
the  forces,  ui  are  the  chemical  potentials,  A ^  are  the  chemi¬ 
cal  affinities  of  the  j  quasi-chemical  reactions,  and  are 
the  reaction  rates.  From  the  second  law  of  thermodynamics, 
the  entropy  source  strength  is  positive  definite  and  equal 
to  zero  in  equilibrium. 

The  basic  assumption  underlying  this  model  is  the  con¬ 
cept  of  "local  equilibrium".  That  is,  although  the  whole 
system  may  be  far  from  equilibrium,  the  particles  in  each 
differential  volume  are  in  equilibrium  with  each  other.  If 
one  considers  the  small  distances  and  interaction  times  among 
particles  within  a  differential  volume  and  the  large  dis¬ 
tances  and  interaction  times  among  particles  from  different 
differential  volumes  of  a  system,  the  concept  of  local  equi¬ 
librium  is  a  reasonable  one. 

For  our  purposes  here  we  neglect  any  temperature  gradi¬ 
ents  or  quasi-chemical  reactions  and  assume  that  the  only 
force  on  the  particles  is  an  electrostatic  potential  gradi¬ 
ent.  Equation  (6)  then  reduces  to, 
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°s  -  -  (7) 

i 

where  is  the  electrochemical  potential.  Consider  the 

form  of  a  near  equilibrium  by  expanding  it  in  a 
s 

series  of  the  generalized  forces  around  =  0.  The  first 

term  is  zero  because  a  *  0  in  equilibrium,  and  the  second 

s 

term,  linear  in  V? . ,  is  zero  because  a  is  positive  definite. 

JL  S 

The  third  term,  quadratic  in  V?.,  is  the  first  non-zero  term. 
Retaining  only  this  term  we  can  see  that,  to  first  order, 
near  equilibrium, 

j±  a-  Vci  .  (8) 


Going  from  particle  current  to  electric  current,  Eq. 
(8)  tells  us  that  the  correct  expression  for  the 
ecruation-of-state  of  this  isothermal  system  is. 


-*■ 


•* 

n. y . V£ . 

ii  i 


(9) 


where  n^  and  y^  are  the  concentration  and  mobility  of  the  i 
particles.  Equation  (9),  of  course,  is  a  well-known  result. 
We  went  through  this  analysis  to  show  that  it  is  correct  to 
first  order,  under  the  assumptions  discussed,  and  that  it 
can  be  derived  from  the  second  law  of  thermodynamics.  In 
general,  all  of  the  quantities  in  Eq.  (9)  can  vary  with 
position  and  time. 

We  now  need  to  obtain  an  expression  for  the  electro¬ 
chemical  potential  for  indistinguishable  particles.  For  a 


grand  canonical  ensemble  of  particles  near  equilibrium,  it 
can  be  shown  tg3  that  the  distribution  function  is 


approximately , 

)f 

That  is,  it  has  the  same  form  as  the  Fermi -Birac  distribution 
except  that,  in  general,  f,  e,  £,  and  T  can  vary  in  time  and 
position.  Equation  (10)  is  obtained  under  the  local  equili¬ 
brium  assumption  and  is  valid  only  to  first  order.  Thus, 

Eqs.  (9)  and  (10)  are  derived  with  the  same  assumptions  and 
approximations,  and  we  can  use  them  to  analyze  heterostruc¬ 
tures  keeping  this  in  mind. 

4.3  HETEROJUNCTION  MODEL 

Let  us  apply  Eq.  (10)  to  obtain  an  expression  between 
the  concentration  of  electrons  in  the  conduction  band  and 
their  electrochemical  potential.  In  the  usual  way,  we 
multiply  the  distribution  function  times  the  density-of-states 
and  integrate  over  the  band  to  obtain, 


g(e)f(e)de  , 


(ID 


where  e  is  the  bottom  of  the  conduction  band.  For  a  dilute 
c 

distribution  of  electrons,  Eq.  (11)  gives. 


n  =  Nc  exp  (C  -  ec) ,  (12) 

where  N  is  the  effective  density-of-states  of  the  conduc- 
c 

tion  band.  In  general,  all  four  quantities  in  Eq.  (12)  can 
vary  with  time  and  position. 


r 
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Using  Eqs.  (9)  and  (12),  we  have  the  general 
near-equilibrium  equation-of-state  for  electron  flow  in  a 
hetero junction, 

3n  “  n“n’(£c  +  kTln  fr)  '  '  <13> 

Notice,  that  all  of  the  quantities  in  this  equation  can,  in 

general,  vary  with  position  and  time  including  N  .  It  is 

c 

the  neglect  of  the  term  -V (  InN  )  which  leads  to  the  incon- 

v 

sistency  mentioned  previously.  A  similar  expression  can  be 
obtained  for  the  hole  current: 


JP  = 


PMpV 


Let's  examine  each  of  the  thermodynamic  forces  in  Eqs. 
(13)  and  (14) .  The  forces  ?ec  and  Vev  represent  variation 
in  composition  or  bandgap  as  well  as  applied  and  built-in 
electric  fields  due  to  space  charge.  The  component  of  these 
forces  which  has  equal  and  opposite  effect  on  electrons  and 
holes  can  be  associated  with  electric  field,  with  the  rest 
associated  with  band-gap  variations.  The  forces  V  (Inn)  and 
-V  llnp)  are  the  usual  diffusion  terms.  The  thermodynamic 
forces  -V ( lhN  )  and  V ( InN  )  can  be  regarded  as  entropy 
forces:  that  is,  they  act  to  drive  the  electrons  and  holes 
to  regions  of  higher  density-of-state  or  higher  entropy. 
These  forces  are  required, therefore ,  by  the  second  law  of 
thermodynamics,  and  must  be  included  in  a  heterostructure 
analysis. 


» 
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We  will  now  show  that  including  these  entropy  forces 
eliminates  the  inconsistancy  in  the  analysis  of  an  abrupt 
p-n  he tero junction.  In  equilibrium,  from  Eqs.  (13)  and  (14) 
the  built  in  potential  of  the  space  charge  region  is. 


Ae  + kTln 
c 


-Aev  +  kTln 


(n  N  \ 
_n_c£\ 

Ncnnp  / 

l  (  PP,Nvn  \ 
\  NvpPn  / 


Using  Eq.  (1)  these  two  expressions  can  be  shown  to  be  equal. 
When  the  junction  is  forward-biased,  the  holes  injected  into 
the  n-side  in  the  depletion  approximation  are. 


Pn  -  Pp  iT*  exp 
*  vp 


q(V-V0)-Aev 


Turning  the  applied  voltage  off  in  Eq.  (17) ,  we  obtain  Eq. 

(16)  and  the  inconsistancy  is  removed. 

Notice  that  in  this  analysis  we  have  made  no  assumptions 
about  the  continuity  of  the  infinite  vacuum  level  (electron 
affinity  rule)  or  any  other  level.  However,  if  we  refer  to 
Eqs.  (7)  and  (8),  in  equilibrium  there  is  no  entropy  produc¬ 
tion  and 


=  o  , 


r 


If  the  electrochemical  potentials  are  referenced  to  the 
infinite  vacuum  level,  then  it  must  be  continuous. 
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5.  INERTIAL  TRANSPORT  DEVICES 

Since  the  primitive  unit  cell  of  the  II-IV-V2  chalco- 
pyrites  is  four  times  as  large  as  that  of  their  III-V 
sphalerite  analogs,  the  first  Brillouin  zone  is  one-fourth 
as  large.  For  this  reason  and  because  of  their  small 
effective  masses,  a  potential  application  for  these 
chalcopyrites  may  be  in  inertial  transport  devices  such 
as  space-charge  limited  diodes  or  Zener  oscillators.  In 
this  section  the  physics  of  these  devices  is  described. 

5.1  SPACE-CHARGE  LIMITED  DIODES 

The  following  reprint  describes  the  current-voltage 


characteristic  of  a  space-charge  limited  diode  in  a 
semiconductor  with  non-parabolic  bands. 
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The  expected  current  is  reduced  significantly  1  -low  the  constant  ef¬ 
fective  mass  value  for  voltages  in  excess  of  the  bjndgap,  and  the  high- 
voltage  asymptotic  relation  is  V  rather  than  V 3I2. 


In  devices  that  use  high-mobility,  low  effective  mass  semi¬ 
conductors,  the  mean  free  time  between  collisions  of  a  carrier 
can  become  equal  to,  or  longer  than  an  oscillation  period,  or 
the  mean  free  path  between  collisions  equal  to  or  longer  than 
a  typical  device  dimension. 

Under  such  conditions  the  carrier  transport  in  semicon¬ 
ductors  is  not  governed  by  the  traditional  friction-dominated 
Ohm’s  law  behavior.  The  motion  of  the  carriers  for  the  short 
times  or  distances  of  primary  interest  for  device  operation  may 
be  more  accurately  described  by  the  inertial  or  ballistic  equa¬ 
tion  of  motion. 

The  simplest  device  for  studying  inertial  transport  is  the 
short  space-charge  limited  semiconductor  diode.  Calculations 
of  the  behavior  of  such  diodes,  without  collisions,  and  with 
few  collisions  [  1  ]  have  been  reported.  These  calculations  as¬ 
sume  that  the  carrier  effective  mass  is  independent  of  energy. 

In  a  device  that  operates  in  a  collision-free,  or  nearly  collision- 
free  regime,  the  carriers  that  contribute  to  the  device  current 
reach  energies  corresponding  to  the  voltages  applied  to  the  de¬ 
vice.  In  small  and  fast  devices,  these  are  of  the  order  of  1  V. 
In  most  semiconductors,  especially  those  with  small  effective 
mass  carriers,  the  bands  depart  from  parabolicity  at  energies 
well  below  1  eV.  The  effective  mass  generally  increases  with 
energy.  Thus  the  current  is  smaller  than  one  would  predict  for 
a  parabolic  band  structure. 

To  assess  the  magnitude  of  this  effect,  we  use  the  conduction 
band  structure  of  a  model  narrow-gap  semiconductor  as  com¬ 
puted  by  two-band  k  •  p  perturbation  theory.  Such  a  model 
can  provide  a  fair  approximation  to  the  real  band  structure  in 
the  vicinity  of  the  conduction  band  extremum. 

Assuming  a  direct-gap  semiconductor,  and  placing  the  zero 
of  energy  at  the  conduction  band  minimum,  the  theory  gives 
the  hyperbolic  band  structure 


1 


2  h2k2 
m*eg 


-  1 


(1) 


Collisionless  Space-Charge  Limited  Currents  in 
Semiconductors  with  Nonparabolic  Bands 


The  theory  also  provides  a  relation  between  the  effective 
mass  m*  and  the  bandgap  eg,  but  for  the  present  purpose  these 
two  parameters  may  be  chosen  independently,  for  instance  to 
fit  a  measured  band  structure. 

To  obtain  the  carrier  velocity  as  a  function  of  energy,  we  use 


9c  d3*  /  |  2 H7k2y/2 
9 k  m*  \  m*eg/ 
/and  then  solve  (1)  for  k 2 


(2) 
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band  structure  as  deduced  from  two-band  k  ■  p  perturbation  theory. 
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Substituting  this  in  (2)  gives 


9c  /~2~  -fl  (e  +  c*/eg)l/2  e/eg  +  e*/e) 
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Fig.  1.  GaAs  band  structure  with  parabolic  and  hyperbolic  conduction 
band  approximation. 


where  we  have  defined  a  band-structure  limited  velocity 

u,  =  (ef/2m*),/J.  (5) 

Fig.  1  shows  the  current  best  plot  of  the  GaAs  band  structure 
(2)  with  both  (1)  using  tne  experimental  6.  and  m*,  and  a 
il^k2  /2m*  parabola  superposed  on  the  plot.  The  hyperbola  of 
(1)  is  evidently  a  much  better  fit. 

The  Poisson  equation  for  a  semiconductor  uniformly  doped 
with  n0  donors/cm3  carrying  an  electron  current  density  J  is 


d2V  p  qn0  ^  J 
dx2  e  e  ev 


(6) 


where  e  is  the  dielectric  constant  and  v  is  the  electron  velocity. 
Electrons  injected  at  x  =  0  acquire  a  velocity 


1  be 
u  =  —  — 
it  bk 


(7) 


when  their  energy  e  -  q  V  s  0. 

Using  the  band  structure  gradient  from  (4)  brings  the  Poisson 
equation  into  the  form 


Jq  (m*u3  +0) 
e  evj  (2 +  <fi2)>/2 


(8) 


where  a  prime  denotes  d/dx.  Equation  (8)  can  be  integrated 
once  to  give 

(0')a  *  (0'(O))J  (<p2  +  2 m*v2s<t>)1'2  -  (9) 

eu,  e 

Assuming  space  charge  limitation  <p\ 0)  =  0,  and  changing  to 
reduced  variables 


_  » 

2m  *v2’ 


qn0vs  ’  (em*ln0q2yl2  v, 


(10) 


we  can  write  the  diode  equation,  the  integral  of  (9),  in  the 
form 


du\i(u2  +  u)1*2  -  u]*1'3 


=  / 


where  u  “  v(V  -  V0,  the  applied  voltage)  at  \-l(x  =  L,  the 


Fig.  2.  Space-charge  limited  current-voltage  relation  of  a  collisionless 
semiconductor. 


diode  length).  The  normalized  variables  are:  diode  voltage  in 
units  of  2m*u|/q— four  times  the  kinetic  energy  the  carriers 
would  have  if  their  effective  mass  were  constant  at  the  band 
edge  value^and  they  were  moving  with  the  limiting  velocity 
if  the  k  ■  p  perturbation  theory  is  taken  literally,  this  unit  is 
equal  to  e./q;  current  density  in  units  of  qnoUj-the  current 
that  would  flow  if  the  mobile  carriers  from  the  doping  were 
moving  with  velocity  v,  \  diode  length  in  units  of  vsTp,  with  Tp 
the  plasma  oscillation  period  of  the  semiconductor. 

Curves  of  current  density  versus  voltage  are  shown  in  Fig.  2 
as  full  curves.  For  comparison  the  current  density-voltage 
curves  one  would  compute  assuming  constant  effective  mass 
are  shown  dashed. 

The  vacuum  diode  analog  3/2  power  law  holds  only  for 
K0  $  e,,  and  then  only  for  li  1.  That  means  very  short 
diodes— for  GaAs  with  n0  =  1016,  shorter  than  0.25  /um.  The 
limiting  higher  voltage  behavior  is  JaV  rather  than  J  a  K3/J; 
the  J-V  curves  bend  over. 

In  the  typical  size  and  operating  voltage  range  of  GaAs  de¬ 
vices  the  corrections  indicated  here  appear  significant  but  not 
dramatic.  They  are  likely  to  be  rather  more  important  in  de¬ 
vices  made  from  smaller  bandgap,  lower  m*  material. 

We  have  not  extended  the  diode  characteristics  into  the  low- 
current  region  at  the  lower  ends  of  the  curves.  The  problems 
encountered  by  modeling  this  operating  regime  are  rather  com¬ 
plicated,  involving  multiple-valued  functions  and  possibly  spa¬ 
tial  and  temporal  instabilities  (3).  Moreover,  the  neglect  of 
collisions  and  thermal  energies  is  less  justifiable  in  this  low- 
voltage  region. 
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5.2  PHASE  FOCUSING  OF  ZENER  OSCILLATIONS 

In  the  absence  of  scattering  and  Interband  tunneling  the  motion 
of  a  carrier  In  a  crystalline  solid  subjected  to  a  steady  uniform  elec¬ 
tric  field  is  described  by  a  rectilinear  trajectory  in  reciprocal  space, 
traversing  the  extended  Brillouin  zone  at  a  constant  rate  proportional 
to  the  field.  The  corresponding  motion  in  direct  space  is  periodic 
with  a  frequency  proportional  to  the  field,  and  is  called  a  Zener  oscil¬ 
lation  T^l  •  Zener  oscillations  are  a  possible  source  of  electromagnetic 
radiation.  They  are  most  likely  to  be  useful  at  extremely  high  frequen¬ 
cies,  in  the  mllllmetrlc  and  submil lime trie  region  of  the  spectrum,  where 
oscillation,  periods  are  comparable  with  or  shorter  than  the  scattering 
lifetime,  so  that  an  appreciable  fraction  of  the  carriers  present  com¬ 
plete  one  or  more  cycles  of  oscillation  before  they  are  scattered. 

The  well-known  equation  of  motion  of  a  carrier  in  a  steady  uniform 
electric,  field  F  is 

hi  -  qZ  (1) 

where  q  is  the  charge  of  the  edrrier  and  hk  its  crystal  momentum.  The 
carrier' s  velocity  v  is 

7  -  i  7tEn(k)  (2) 

where  Eq  is  the  momentum-dependent  energy  in  the  band  of  index  n.  If 
k  lies  along  a  principal  crystal  direction,  the  period  of'  the  Zener 
oscillation  is 


T 


t 


(3) 


where  <  is  the  diameter  of  the  first  Brillouin  zone  in  the  direction 
k,  and  the  length  of  the  Zener  orbit  is 
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where  E  _  and  E _ are  respectively  the  maximum  and  minimum  energy 

nmax  ntnm 

attained  by  the  carrier  in  the  course  of  its  traversal  of  the  Brlllouin 
zone. 

For  a  typical  cubic  semiconductor  of  lattice  constant  of  order  5$, 

in  an  electric  field  of  order  10^V/cm  the  Zener  oscillation  frequency  is 
12 

of  order  10  Hz,  and  if  the  energy  width  of  the  band  is,  say,  SeV,  the 
length  of  the  Zener  orbit  is  0.5  pm.  (It  is  convenient  to  remember  that 
the  length  of  the  Zener  orbit  is  equal  to  the  distance  in  which  the 
electrostatic  potential  of  the  applied  field  drops  by  the  energy  width 
of  the  band.) 

These  equations,  numbers,  and  statements  are  to  be  understood  in 
terms  of  a  quasi -classical  picture  of  the  crystal  electrons.  This  picture 
is  commonly  accepted  for  the  heuristic  description  of  carriers  in  solid 
state  electron  devices,  and  with  appropriate  caution  it  can  be  used 
here.  It  is  well  established  that  interband  tunneling  is  negligible 
in  applied  potentials  which  differ  by  much  less  chan  the  bandgap  at 
adjacent  lattice  sites,  in  the  sense  that  the  tunneling  probability  per 
Zener  period  is  much  less  than  unity.  We  shall  show  elsewhere  that  in 
such  fields  crystal  electrons  can  indeed  be  localized  simultaneously  in 
momentum  in  a  region  s^iall  compared  with  a  Brlllouin  zone,  and  in  posi¬ 
tion  in  a  region  small  compared  with  the  Zener  orbit  [5] .  The  localiza¬ 
tion  is  not  Independent  of  time,  but  like  a  coherent  state  of  a  har¬ 
monic  oscillator  [6]  (and  unlike  a  wave  packet  representing  a  free  elec¬ 
tron)  it  does  not  spread  indefinitely  but  recurs  with  the  Zener  period. 
This  property  arises  from  the  equal  energy  spacing  of  the  Stark  ladder 
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states  that  constitute  the  wave  packet.  We  also  accept  here  as  plaus¬ 
ible  (and  shall  demonstrate  elsewhere)  that  a  small  change  In  the  field 
orientation  only  slightly  affects  the  Zener  spectrum  [5]. 

**  The  power  radiated  by  oscillating  charges  Is  small  unless  they 

oscillate  In  phase.  A  striking  If  trivial  example  Is  presented  by  the 
'  charges  In  filled  bands  of  a  solid.  If  a  field  is  applied  to  the  solid, 

each  of  these  charges  moves  according  to  Eq.  (1) ,  executing  Zener  os¬ 
cillations  (without  being  scattered,  because  all  the  final  states  are 
full).  Of  course,  these  oscillations  are  electromagnetlcally  unobserv- 
able,  since  the  "radiated"  fields  of  all  the  charges  cancel  because  of 
their  uniform  distribution  over  2ir  radians  of  phase. 

In  order  to  use  the  Zener  oscillations  of  carriers  in  a  partially 
P  filled  band  for  the  coherent  generation  of  radiation.  It  Is  necessary 

to  start  all  the  carriers  oscillating  together  -  phase  Initialization  - 
and  to  maintain  them  in  step  as  they  oscillate  -  phase  focusing. 

We  propose  to  meet  both  these  requirements  by  means  of  a  device 
which  is  schematically  represented  in  Fig.  1.  In  Fig.  1(a)  is  shown 

/ 

band  diagram  typical  of  a  semiconductor  heterojunction  quantum  well. 

Such  planar  structures  have  been  fabricated  by  heteroepitaxy  in  a 
variety  of  semiconductors  [7].  They  can  be  designed  by  varying  the 
relative  stoichiometry  and  doping  of  the  constituent  materials  and  the 
—  thickness  of  the  layer  comprising  the  well  to  contain  one  or  a  few 

one-dimensional ly  bound  electron  and/or  hole  states,  as  indicated  in 
the  figure.  In  Fig.l(b)the  band  diagram  of  the  same  structure  is 
9  shown  with  a  uniform  steady  electric  field  applied  at  right  angles  to 

the  plane  of  the  epitaxial  layer.  In  the  presence  of  such  a  field, 
it  becomes  energetically  possible  for  the  bound  electron  to  tunnel 
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Conduction  Band 


Valence  Band 


Figure  A.  Quantum  Well  Structure 


Figure  B.  Quantum  Well  In  Transverse 
Electric  Reid. 


Figure  1  Schematic  band  diagram  of  a  heteroj unction  quantum  well: 
(a)  in  equilibrium;  (b)  with  an  applied  transverse 
electric  field. 


through  the  energy  barrier  confining  it  to  the  quantum  well,  i.e.  to 
be  Internally  field-emitted  into  the  conduction  band  of  the  "host" 
semiconductor,  to  the  right  of  the  quantum  well  in  the  figure.  Once 
the  electron  has  been  field-emitted  into  the  conduction  band,  it  will 
execute  Zener  oscillations  until  it  encounters  a  scattering  event. 

Upon  being  scattered,  the  electron  will  resume  its  oscillatory 
motion,  in 'general  with  a  new  phase  unrelated  to  the  phase  prior  to 
the  scattering  event,  thus  becoming  useless  or  harmful.  But  observe 
that  the  motion  of  these  dephased  electrons  in  general  resumes  "down¬ 
hill"  from  the  original  starting  point,  so  that  their  Zener  oscillation 
orbits  extend  to  the  right  beyond  the  maximum  excursion  of  the  favorably 
phased  carriers.  If  an  appropriate  collecting  structure  is  placed  in 
this  region  a  distance  just  in  excess  of  l  in  Eq.  (4)  away  from  the 
quantum  well  it  will  eliminate  the  dephased  carriers  from  the  inter¬ 
action  region.  This  constitutes  the  phase  focusing  mechanism  of  the 
device  concept.  From  the  parenthetical  remark  made  above  in  connection 
with  Eq.  (4)  it  follows  that  the  length  of  the  Interaction  region  be¬ 
comes  equal  to  t  when  the  voltage  across  it  reaches  q  ^(EnniaY  -  Enni1n) ; 
thus  the  spacing  of  the  collecting  electrode  sets  the  lower  frequency 
limit  of  the  device  through  Eqs.  (4)  and  (3). 

It  is  convenient  to  discuss  the  oscillation  of  the  planar  charge 
density  and  the  associated  time-varying  electromagnetic  fields  in 
lumped-circuit  terms.  In  this  language,  the  continuity  of  the  total 
(conduction  and  displacement)  current  associated  with  the  oscillating 
charge  requires  a  current  to  flow  in  any  external  circuit  surrounding 
it.  This  current  in  turn  will  induce  a  voltage  that  alters  the  field 
in  the  interaction  region.  If  the  external  circit  Impedance  is  a 


capacitance,  the  phase  of  this  voltage  is  such  as  to  increase  the  applied 
steady  voltage  at  those  times  when  we  wish  to  have  carriers  field-emitted 
across  the  potential  barrier  forming  the  edge  of  the  quantum  well  in  order 
to  oscillate  in  synchronism  with  carriers  field-emitted  during  earlier 
cycles  of  oscillation.  Since  the  tunneling  probability  depends  very 
sensitively  on  the  width  of  the  potential  barrier,  and  since  the  width 
is  smallest  when  the  voltage  is  largest,  the  feedback  from  the  external 
circuit  serves  to  maintain  coherence  in  the  phase  Initialization. 

The  capacitive  load  required  for  the  appropriate  phasing  of  the 
feedback  can  be  supplied  in  a  variety  of  ways.  Because  of  the  high  * 
operating  frequency  of  the  device,  it  may  be  convenient  to  use  quasi¬ 
op  t  leal  circuitry  together  with  a  junction  capacitance.  Alternatively, 
a  microwave-type  implementation  could  have  the  device  operate  into  a 
short-circuited  transmission  line  between  a  quarter  and  a  half  wavelength 
long. 

The  collecting  structure  should  be  an  interface  with  an  effectively 
infinite  recombination  velocity.  This  suggests  a  reverse-biased 
semiconductor-metal  junction  as  the  positive  electrode.  Which  also  serves 
as  the  bias  supply. 

Fig.  1(b)  has  been  drawn  with  the  assumption  that  the  field  in  the 
Interaction  region  is  uniform.  This  Implies  that  the  semiconductor  is 
intrinsic.  Design  considerations  may  make  it  advantageous  to  use  doped 
material.  In  this  case  the  field  in  the  interaction  region  would  be 
non-uniform,  implying  trivial  modification  of  the  equations.  Whether 
intrinsic  or  doped,  the  interaction  region  is  depleted  under  operating 
conditions,  that  is  to  say,  the  device  is  operated  in  punch-through. 


Another  simplifying  feature  of  Fig.  1  is  the  symmetry  of  the  walls 
of  the  quantum  well,  by  no  means  a  necessary  condition.  By  varying  the 
doping  on  the  "upstream"  side.  It  Is  possible  to  control  the  design 
and  the  high-frequency  circuit  function  of  the  junction  at  the  nega¬ 
tive  contact  which  resupplies  carriers  to  the  well. 

In  the  mode  described  here  the  proposed  device  will  function  best 

at  frequencies  above  the  scattering  rate,  which  Is  typically  of  order 
12 

10  /sec,  that  Is  to  say.  In  the  Terahertz  region  of  the  spectrum.  As 

seen  In  the  numerical  estimate  following  Eq.  (4),  the  dimensions  and 

12 

operating  parameters  of  a  10  Hz  device  are  quite  reasonable.  Operation 
at  lower  frequencies  should  be  possible,  since  the  proposed  phase-focusing 
mechanism  eliminates  scattered  carriers.  The  efficiency,  however,  would 
suffer,  since  the  scattered  carriers  contribute  to  the  total  current  and 
hence  to  the  Input  power. 

The  scattering,  of  course,  also  contributes  to  the  heat  Input  to 
the  device.  One  may  expect  a  lower-frequency  limit  to  the  operation  to 
be  set  by  this  heating.  Since  a  lower-frequency  device  also  needs  to 
be  thicker.  It  may  reach  the  avalanching  threshold,  also  setting  a 
low-frequency  limit. 
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5.3  TERAHERTZ  ZENER  OSCILLATOR 
5.3.1  Device  Concept 

Zener  oscillations  are  the  real-space  oscillations  that  correspond 
to  the  uniform  rectilinear  reciprocal-space  trajectory  of  a  crystal  elec¬ 
tron  in  a  steady  applied  electric  field.  Such  oscillations  have  a  period 
r  given  by 

r  ■  ffic/qE 
and  an  amplitude  t  given  by 

t  -  (1/qE)  -  V„) 

where  q  is  the  electronic  charge,  E  the  electric  field,  <  the  "diameter” 

of  the  Brillouin  zone  along  the  field  direction,  and  E  and  E  ,  respect- 

mix  min 

ively  the  maximum  and  minimum  of  the  band  structure  traversed  along  the 
electron's  trajectory.  These  expressions  assume  the  absence  of  scatter¬ 
ing  and  interband  tunneling  events  during  time  intervals  of  interest,  as¬ 
sumptions  to  be  discussed  below. 

Oscillations  of  this  type  are  capable  of  generating  high-frequency 
electromagnetic  radiation.  Zf  the  charge  carriers  are  to  radiate  effici¬ 
ently,  they  should  oscillate  in  phase.  We  propose  a  scheme  whereby  car¬ 
riers  in  a  semiconductor  are  injected  into  a  depleted  interaction  region 
Periodically,  with  the  period  of  the  Zener  oscillations,  and  are  automat¬ 
ically  removed  if  they  are  scattered  and  lose  phase  coherence. 

The  principle  of  the  proposed  device  is  schematically  illustrated 
in  Fig.  1.  Fig.  la  is  the  thermal  equilibrium  band  diagram  of  a  structure 
consisting  of  a  wide-bandgap  intrinsic  semiconductor  interaction  region 
between  a  narrower-gap  n-doped  semiconductor  on  the  left  and  a  metal  con¬ 
tact  on  the  right.  In  Fig.  lb  the  structure  is  shown  with  a  positive 
voltage  applied  to  the  metal.  It  is  now  possible  for  electrons  to  tunnel 
through  the  n-i  hetero junction  barrier  (to  be  internally  field-em  .tted) 
into  the  interaction  region,  where  they  will  execute  Zener  oscillations. 


A  moment's  thought  reveals  that  the  length  of  the  Zener  orbit  H  becomes 

equal  to  the  width  of  the  interaction  region  when  the  voltage  across  it 

is  equal  to  Vq  *  (1/q)  (smjty  ”  n)  .  Thus  with  a  voltage  equal  to  or 

greater  than  V  across  the  interaction  region,  a  conduction  electron's 
c 

Zener  orbit  remains  confined  to  that  region,  unless  it  experiences  a 
scattering  event. 

If  it  is  scattered,  the  electron  will  resume  its  oscillatory  notion, 
in  general  with  a  new  phase  unrelated  to  the  phase  prior  to  the  scattering 
event,  thus  becoming  useless  or  harmful.  But  observe  that  the  motion 
of  these  dephased  electrons  in  general  resumes  "downhill"  from  the  ori¬ 
ginal  starting  point,  so  that  their  Zener  orbits  extend  to  the  right, 
beyond  the  maximum  excursion  of  the  favorably  phased  carriers.  Thus  if 

the  voltage  across  the  interaction  regisn  exceeds  v  only  slightly,  these 

c 

dephased  carriers  will  be  collected  at  the  metal  contact  and  eliminated 
from  the  interaction  region. 

The  oscillating  planar  charge  density  of  conduction  electrons  and 
the  associated  time-varying  electric  fields  constitute  a  current.  The 
total  {convection  plus  displacement)  current  must  be  continuous  across 
the  junctions  and  in  the  external  re turn  circuit.  The  current  in  the 
external  circuit  in  turn  induces  a  voltage  which  appears  across  the  inter¬ 
action  region.  If  the  reactance  of  the  external  circuit  is  a  capacitance, 
the  phase  of  this  voltage  is  such  as  to  increase  the  applied  steady  volt¬ 
age  at  those  times  when  we  waext  the  n-region  to  inject  new  carriers  into 
the  interaction  region  to  oscillate  in  synchronism  with  carriers  emitted 
during  earlier  cycles.  Since  the  tunneling  probability  depends  very  sens¬ 
itively  on  the  width  of  the  hetero junction  potential  barrier,  and  since 
thewidth  is  smallest  tAien  the  voltage  is  largest,  the  feedback  from  the 
external  circuit  serves  to  maintain  coherence  in  the  phase  initialization. 

The  feedback  loop  is  analogous,  in  some  measure,  to  that  of  an  IMP ATT 
oscillator.  The  transit-time  phase  delay  is  similar  in  both  devices.  The 
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IMP  ATT'  s  delay  due  to  the  avalanche  multiplication  is  replaced  here  by 
the  phase  shift  of  the  circuit  reactance. 

5.3.2  Discussion 

The  description  of  the  proposed  device  has  been  couched  in  terms  of 
a  quasi-classical  picture  of  the  crystal  electrons.  This  picture  is  com¬ 
monly  accepted  for  the  heuristic  description  of  carriers  in  solid-state 
electron  devices,  and  with  appropriate  caution  it  can  be  used  here.  It 
is  well  established  that  interband  tunneling  is  negligible  in  applied 
potentials  which  differ  by  much  less  than  the  bandgap  at  adjacent  lattice 
sites,  in  the  sense  that  the  tunneling  probability  per  Zener  period  is 
much  less  them  unity. 

We  have  also  assumed  that  the  crystal  electron  can  be  localized  sim¬ 
ultaneously  in  momentum  in  a  region  small  compared  with  the  Brillouin 
zone,  and  in  position  in  a  region  small  compared  with  a  Zener  orbit. 

Such  localization  is  routinely  taken  for  granted  in  low  fields,  but  must 
be  examined  more  carefully  in  the  present  context.  For  reasons  to  be 

discussed  presently,  we  shall  assume  a  Zener  oscillation  frequency  in 
12 

excess  of  10  Hz.  In  a  typical  cubic  semiconductor  of  lattice  constant 

o  5 

of  order  5A,  this  requires  a  field  of  order  10  V/cm,  and  if  the  energy 
width  of  the  band  is,  say,  5eV,  the  length  of  the  Zener  orbit  is  0.5ym. 
we  have  constructed  gaussian  wavepackets  of  one-band  solutions  of  the 
Schroedinger  equation  (Stark  and  Houston  states)  for  a  model  tight-bind¬ 
ing  band  structure,  and  verified  that  in  a  field  of  10^V/cm  such  wave 
packets  can  be  simultaneously  confined  to  less  than  1%  of  the  Brillouin 
zone  and  to  less  than  1/30  of  a  Zener  orbit  length.  The  localization  in 
position  is  not  independent  of  time,  but  like  a  coherent  state  of  a  har¬ 
monic  oscillator  f  9]  (and  unlike  a  wave  packet  representing  a  free  electron) 
it  does  not  spread  indefinitely  but  recurs  with  the  Zener  period.  This 
property  arises  from  the  equal  energy  spacing  of  the  Stark  states  that 


constitute  the  wave  packet.  (In  the  numerical  example  given  above,  the 
position  uncertainty  fluctuates  by  about  50%  in  the  course  of  a  cycle  of 
oscillation.) 

Some  concern  has  been  expressed  in  the  literature  about  the  stability 
of  the  Stark  ladder  spectrum  "in  the  sense  that  the  slightest  change  in 
the  direction  of  E  completely  alters  the  level  system"  [9].  Although  this 
statement  is  strictly  true  -  indeed,  the  Zener  oscillations  are  aperiodic 
unless  E  lies  along  a  reciprocal  lattice  vector  -  it  has  little  bearing 
on  our  present  concerns.  Using  a  realistic  GaAs  conduction  band  structure  [10] 
we  have  determined  that  the  effect  of  a  small  change  in  the  direction 

of  E  is  akin  to  a  slow  frequency  modulation  of  the  Stark  spectrum. 

12 

.  We  have  chosen  a  Zener  frequency  of  order  10  Hz  and  a  correspond¬ 
ing  field  of  order  l05V/cm  as  appropriate  for  the  testing  of  our  device 
concept.  This  choice  is  supported  by  the  experimental  results  of  Koss 
and  Lambert  [11]  ,«ho,in. fields  of  this  order  of  magnitude,  claimed  to  ob¬ 
serve  the  Stark  ladder  in  the  GaAs  conduction  band.  The  result  consis¬ 
ted  of  a  step-like  modulation  of  the  Franz-Keldysh  effect  on  the  interband 
optical  absorption.  Although  the  resolution  of  the  experiment  was  not 
very  good,  and  thus  the  evidence  is  not  entirely  conclusive,  it 
implies  that  the  lifetime  of  a  Stark  level  is  of  order  10-1,2s. 

A  lifetime  of  this  order  of  magnitude  appears  to  be  in  conflict  with 
a  variety  of  experimental  and  theoretical  results  [12]  on  high-field  trans¬ 
port  and  avalanche  ionization  which  yield  scattering  rates,  mainly  assoc¬ 
iated  with  phonon  emission,  of  at  least  10^*/s.  Although  we  are  unable 
«t  present  to  account  for  this  apparent  discrepancy,  we  can  suggest  a 
direction  in  which  a  resolution  might  be  found. 

Existing  high-field  transport  theory  is  based  on  calculations  involv¬ 
ing  the  interaction  of  quasi-particles  represented  as  plane  waves.  Such 
theories  are  strictly  valid  only  in  effectively  infinite  homogeneous  media. 
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Although  for  many  experimental  situations  this  is  a  quite  reasonable  ap¬ 
proximation,  it  may  be  expected  to  fail  when  the  scale  of  the  non-uniform¬ 
ities  becomes  small  enough.  In  particular,  in  the  device  proposed  here 
we  deal  with  very  tightly  localized  quasi-particles.  This  localization 
must  be  at  least  partially  conserved  in  a  scattering  event.  One  may  ex¬ 
pect  that  this  requirement  places  a  constraint,  in  effect  a  set  of  select¬ 
ion  rules,  on  the  type  of  scattering  event  that  may  occur.  As  is  gener¬ 
ally  true  of  any  constraint,  the  restrictions  thus  imposed  must  reduce 
the  probability  of  occurrence  of  the  event  and  increase  the  lifetime. 
Moreover,  the  allowed  events  may  have  less  drastic  experimental  conse¬ 
quences  than  a  complete  randomization.  For  example,  coherent  forward 
scattering,  even  with  the  loss  of  a  phonon  energy,  may  not  be  fatal  to 
the  usefulness  in  the  interaction  region  of  a  Zener  oscillating  carrier. 
Such  qualitative  arguments  have  a  good  deal  of  plausibility,  to  render  them 
quantitative  will  require  methods  that  have  recently  been  developed  to  deal 
with  high-field,  small-device  transport  phenomena  [13  J. 

Both  the  magnitude  and  the  field  sensitivity  of  the  hetero junction 
tunneling  current  will  be  satisfactory  for  this  application.  These  quant¬ 
ities  can  be  estimated  using  the  Fowler— Nordheim  equation.  With  a  barrier 
height  of  a  few  hundred  meV  and  a  field  of  10  V/cm  we  find  current  densi¬ 
ties  comparable  with  those  of  tunnel  diodes.  The  fractional  field  sen¬ 
sitivity  is 


E  dl 


I  dF  "  2  +  3ft  qF 
where  I  is  the  tunneling  current,  m  the  carrier  effective  mass,  and  H  the 
barrier  height.  This  ratio  can  be  adjusted  by  varying  the  doping,  and  it 
will  generally  exceed  10. 

The  capacitive  load  required  for  the  appropriate  phasing  of  the  feed¬ 
back  can  be  designed  in  a  variety  of  ways.  Because  of  the  high  operating 
frequency  of  the  device,  it  may  be  convenient  to  use  quasi-optical  cir- 


4/(2mH)  H 


cuitry .  Then  the  capacitance  might  be  provided  by  a  reverse-biased  junct¬ 
ion  on  the  narrower-gap  side  of  the  hetero junction.  Alternatively,  a 
microwave- type  implementation  could  have  the  device  operate  into  a  short- 
circuited  transmission  line  between  a  quarter  and  a  half  wavelength  longr. 

Fig.l  has  been  drawn  with  the  assumption  that  the  interaction  region 
is  intrinsic.  This  is  not  a  necessary  condition,  and  design  consider¬ 
ations  may  make  it  advantageous  to  use  doped  material.  Zn  this  case  the 
field  in  the  interaction  region  would  be  non-uniform,  implying  trivial , 
modification  of  the  equations.  Whether  intrinsic  or  doped,  the  inter¬ 
action  region  under  operating  conditions  is  depleted,  that  is  to  say,  the 
device  is  operated  in  punch- through. 
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